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Recently, three-dimensional Stirling engine simulations have been accomplished utiliz- 
ing commercial Computational Fluid Dynamics software. The validations reported can be 
somewhat inconclusive due to the lack of precise time accurate experimental results from 
engines, export control/proprietary concerns, and the lack of variation in the methods 
utilized. The last issue may be addressed by solving the same flow problem with alter- 
nate methods. In this work, a comprehensive examination of the methods utilized in the 
commercial codes is compared with more recently developed high-order methods. Specif- 
ically, Lele’s Compact scheme and Dyson’s Ultra Hi-Fi method will be compared with 
the SIMPLE and PISO methods currently employed in CFD-ACE, FLUENT, CFX, and 
STAR-CD (all commercial codes which can in theory solve a three-dimensional Stirling 
model although sliding interfaces and their moving grids limit the effective time accuracy). 

We will initially look at one-dimensional flows since the current standard practice is to 
design and optimize Stirling engines with empirically corrected friction and heat transfer 
coefficients in an overall one-dimensional model. This comparison provides an idea of the 
range in which commercial CFD software for modeling Stirling engines may be expected 
to provide accurate results. In addition, this work provides a framework for improving 
current one-dimensional analysis codes. 

Nomenclature 


/3 Nondimensional Wavenumber 

e Turbulent Dissipation Rate, (?n 2 /s 3 ) 

p Kolmogorov Length Scale, (m) 

/i Dynamic (Molecular) Viscosity, (As • m~ 2 ) 

v Kinematic Viscosity, (p • p _1 ) 

p Density, (kg ■ m~ 3 ) 

t w Surface Shear Stress, (kg ■ m~ 1 s~ 2 ) 

c Constant Convective Velocity 

c4o0 UHF Method - 4 point stencil - no derivatives 

c4ol UHF Method - 4 point stencil - one derivative 

c4o2 UHF Method - 4 point stencil - two derivatives 

c4o3 UHF Method - 4 point stencil - three derivatives 

C p Specific Heat Constant Pressure 

C v Specific Heat Constant Volume 

G Complex Amplification Factor 

k Wavenumber 

L Length of Numerical Domain 

r Courant Number 

u T Friction Velocity, (kg ■ m~ 3 ) 

v Von Neuman Number 

CFL Courant-Friedrichs-Lewy 
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I. Introduction 


P OWER conversion with free-piston Stirling engines 1 promises to deliver high efficiency, low mass solutions 
for longer and more varied space missions. 2 In addition to using advanced high-temperature materials to 
increase the Carnot temperature ratio, it is anticipated that advanced computational fluid dynamics (CFD) 
will help to identify the following losses 4, 5,6 (also shown in figure 1): 

1. Inefficient heat exchange and pressure loss in the heat exchangers (heater, regenerator, and cooler) 

2. Gas spring and working space loss due to hysterisis and turbulence, 

3. Appendix gap losses due to pumping and shuttle effects, 

4. Mixing gas losses from unequal temperature distributions or losses from mixing of gas streams, or 
elements of gas at different temperatures, 

5. Conduction losses from the hot to cold regions 

In addition, the following artificial numerical losses must be considered when computational simulations 
are performed (also shown in figure 1): 

1. Moving/ deforming mesh losses from repeated low order flow field interpolations, 

2. Transient/Unsteady heat transfer and flow loss from inconsistent and inaccurate time discretization, 

3. Flow loss from low order approaches resulting in effectively adding artificial dissipation terms along 
sliding interfaces, at structured/unstructured grid interfaces, and within interior. 

Minimizing those artificial losses is best accomplished through higher order approaches.' While this 
approach is common in aeroacoustics, computational electromagnetics, and exterior flow problems, high 
order techniques have not yet been applied to simulating a Stirling device. Moreover, the following difficulties 
are often encountered when using high-order approaches: 

1. Generation of high-order, smooth, body-fitted grids around complex configurations can be difficult. 8 

2. High-order formulations can lack nonlinear robustness. 

3. The general usefulness of high-order methods is limited by first order accurate shock capturing. 9 

Fortunately, with the exception of the possibly random geometry in the regenerator, the free-piston design 
is essentially smooth and admits curvilinear structured grids (with some geometry simplification) . The issue 
of nonlinear robustness (i.e. maintaining design accuracy with nonlinear equations) is an open issue, but 
preliminary results are encouraging. And finally, the working gas is subsonic and shockless throughout the 
entire region 10 (However, steep temperature gradients can exist at solid/fluid interfaces). For these reasons, 
a high-order approach is investigated for ’’whole engine” Stirling analysis and compared to commonly used 
techniques. 


II. Description of the Problem 

The dual opposed configuration shown in figure 2 11,12 is being developed for multimission uses. 13 

Many methods in general use stop at A th order accuracy for time dependent problems since they use 
Runge-Kutta methods. High-order Runge-Kutta methods become notoriously difficult to derive because the 
number of nonlinear order conditions that need to be solved grows exponentially (i.e., a 12 th order method 
has 7813 nonlinear order conditions). The advantages of using Runge-Kutta methods at orders less than 6 
are commonly cited as flexibility, large stability limits, and ease of programming. 14 The practical limit on 
their order has been an impediment to the analysis of their use in high order approaches for time dependent 
applications. 15 

In this paper we use a series of explicit, local, high order methods which have the same order of accuracy 
in space as in time 16,17 for inviscid flow (lower accuracy in time for viscous flows). These methods use 
Hermite interpolation on stencils that are four points wide, and a Cauchy-Kowalewski recursive procedure 18 
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Schematic with Springs and Dampers 


Effect of Heating 



Bearing System 


Moving Mesh Loss 


& Cooling 



Artificial Transient Heat Loss 


Flow Losses (or Artificial CFD Losses) 


Gas Spring 


Piston/Displacer Coupling via Flow Losses 

X„ 


Load 

Working Fluid (bounce space and working space) 


Figure 1. Schematic of Actual and Artificial Numerical Losses in a Free-Piston Convertor 


for obtaining time derivatives from the space derivatives of the interpolant. The time derivatives are then 
used to advance the primitive variables and their spatial derivatives in time with a Taylor series expansion. 
This general approach is called the Modified Expansion Solution Approximation (MESA) method 19 and the 
new finite volume variation of this is called the Ultra HI-FI (UHF) method. 20 This method can be used 
to derive and implement algorithms with arbitrarily high orders of accuracy in multiple space dimensions if 
their complexity is properly managed and the computer’s floating point precision is sufficiently high. 21 

First, some of the known exact solutions of the viscous Burger’s equation are provided and the linear 
case is solved with both state-of-practice Compact schemes, new UHF methods, and various commercial 
code solvers. The one-dimensional Navier-Stokes equations reduce to the linear viscous Burger’s equation 
in certain circumstances and provide a means for testing both heat transfer effectiveness and turbulent 
transition efficiency of each method. 



Figure 2. Dual Opposed Stirling Convertors Reduce Vibration 
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III. Exact Solutions For Method Comparison 

Since the nonlinear Navier-Stokes equations generally clo not have exact solutions nor known stability 
limits, preliminary development and testing of new numerical methods is best accomplished by starting with 
the viscous Burger’s equation. This equation describes flow behaviour in specialized circumstances, but more 
importantly, its mathematical properties are very similar to the full Navier-Stokes equations and it admits 
exact solutions. 

For future reference and convenience, some of the known exact solutions are shown below (only the linear 
viscous Burger’s equation will be required in this work). 

A. Complete Nonlinear Viscous Burger’s Equation 

The viscous Burger’s equation is written as: 


dt U dx ^ dx 2 

where u is the convective velocity term and /x can be considered the dynamic viscosity. 
Exact steady-state solution, hm t ^ (X) u(a;, t) exists for the case with boundary conditions: 


(1) 


1l(0, t) = Uq 
u(L, t) = 0 


and it is given by: 


1 — exp[uReL(x/L — 1 )] 
0 1 + exp[ui?ez,(a:/L — 1)] 

where Re^ = (note this is a modified Reynold’s number) 
u is a solution of the equation 

— = exp (~uRe L ) 
u + 1 

Other solutions for the nonlinear viscous burgers equation are: 


1 . 


u t + uu x - fm xx = 0, /x = 0.1, x G (0, 1) 


(2) 

( 3 ) 

( 4 ) 

( 5 ) 

( 6 ) 


u(a:,0) =0,w(l,t) = — tank 



, u( 0, t) = 0 


has exact solution 


u(x, t) = — tank 



2 . 


Fully nonlinear equation: 


u t + uu x - ~{uu x ) x 

with initial and boundary conditions as: 

u(x, 0) = exp^, u(0, t) 

has the following exact solution u(x,t) = exp®. 


0,x G (0,1) 
l,it(l,f) = e 


( 7 ) 


(8) 

(9) 


(10) 
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3. An exact solution of the nondimensional form of the Nonlinear Burger’s equationn 22 


Nondimensionalized by: 


producing: 


One stationary solution is: 


du du d 2 u 
= »- 


dt 


dx 


dx 2 


* x * uL n 

x = —,u = — ,t* = 

L n L z 


u* * du* d 2 u* 

- 1 - 'll = 

t* dx* dx* 2 


2 sinh x* 


u = — 


cosh x* — exp 


( 11 ) 


(12) 


(13) 

(14) 


B. Linear Viscous Burger’s Equation 

For simplicity and more thorough stability analyses, the viscous Burger’s equation may be linearized with 
constant convective velocity, c, and dynamic viscosity, /i: 


du 

dt 


du 


d 2 u 
dx ^ dx 2 


(15) 


The exact steady-state solution with the same boundary conditions as in Eqs.(2) and (3) is: 


u = Uq 


1 - exp [R l (x/L - 1)] 
1 - exp(-i? L ) 


where Rl = — (modified Reynolds number). 

The exact solution for the linearized equation with initial condition, u(x, 0) 
boundary conditions is: 

u(x, t) = exp (— fc 2 /it) sin k(x — ct) 


(16) 


sin(kx), and periodic 


(17) 


This is useful for evaluating the temporal accuracy of a method and this will be used in comparing Compact 
and UHF techniques . 23 

In addition, this equation form also describes the time-accurate temperature distribution in a moving 
solid or within a moving fluid in a channel in which case fi = a = is interpreted as the thermal diffusivity, 
and T is the temperature: 


Q T g T d 2 T 

~dt + U ~dx ~ a ~dx 2 ^ ^ 

This will be used for comparing the heat transfer capabilities of Compact, UHF, and commercial code 
solvers. 


C. 2-D Nonlinear Viscous Burgers’ Equation 

The 2-D Burger’s equation with nonlinear convection terms , 24 


with 


du du du ,d 2 u d 2 u 

dt + U frc + V &y ~ + fry 2 ) ~ fx = 


dv dv dv ,d 2 v d 2 v. 

u — + v — - + lOi) “ fy = 0 


dt 


dx 


dy 


' dx 2 dy 2 ' 


„ 1 x 2 + 2xy „ o 9 

fx = — 7 H — t + 3 x 3 y — 2 vy 

Jx (1 + t) 2 (1 + t) y y 


(19) 

(20) 


(21) 
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+ 3 y 3 x 2 — 2vx 


(22) 


1 y 2 + 2 xy 
(1 + t ) 2 (1 + t) 

has exact solution: 

u = + x 2 y (23) 

v = + x v 2 ( 24 ) 


D. 2-D Linear Burger’s Equation 

The 2D (Burger’s) linear convection and diffusion equation: 25 

■Ut + c(u x + Uy) - n{u xx + Uyy) = 0, (x , y) G (-1, 1) x (-1, 1); c = 1, y = 0.01 (25) 

with the initial condition u(x, y , 0) = sin(7r(a;+j/)) and periodic boundary condition has the exact solution 
is: 

u(x, y , t) = exp _2lrAlt sin(7r(a; + y — 2 cf)) (26) 

Additional exact solutions may be found in the paper by Benton and Platzman. 26 


IV. Application of Compact Scheme 


The currently accepted state-of-the-art approach to high fidelity numerical simulations is based on Lele’s 
6 th order Compact scheme. 27 This approach implicitly solves the spatial derivative terms and utilizes stan- 
dard Runge-Kutta time advance. Each time step requires solving a tridiagonal matrix given by the equation: 


a 



f du\ ( du\ _ U i+ i - Ui-i U i+ 2 - Ui-2 

\dx ) t a \dx) i+1 a 2 Ax 4Ax 


with a = 1/3, a = 14/9, b = 1/9. 

Similarly, the second order derivatives are found by: 


(27) 


a 



( d 2 u\ ( d 2 u\ _ u i+1 - 2m + Uj-i u i+ 2 - 2u.j + Uj -2 

\dx 2 ) i a \dx 2 ) i+1 a Ax 2 4Ax 2 


with a = 2/11, a = 12/11, b = 3/11. 

The 4 th order Runge-Kutta method is described by: 24 


(28) 


R(u) 

— Cl lx H - l^xx 

(29) 

U V> 

= u n + ^-R n 

(30) 

u™ 

= 

(31) 


= u n + AtR 2 

(32) 

u n+1 

= u" + ^ ( R n + 2 i? (1) + 2 i? (2) + f? (3) ) 

(33) 


(34) 


with 

f? (1) = R( u ^) 
R {2) = f?(u (2) ) 
f? (3) = R(u (3) ) 
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The single step 6 th order Compact scheme with 4 th order Runge-Kutta on a domain [—4, 4] is given by: 


n+1 


( 6296179542799261293 r 4 
^ 2934584609881498112 


98816474558753811m- 3 

18549839506204160 


3 1068 146188 125r 3 
42158726150464 


3588625 ur 


287760797433ot 2 

133245025760 

573327»r 


3898048 


1093840 


75r 

8701 


324338625r 2 , 401096399245idr 
2422636832 “ l " 

_ 1230500161v 4 
24196548096 


583491009024 
268897'0 3 , 16675t> 


55178557611^ r 
19381094656 

It— 4 T 


It— 4 + 


1714608 


127008 


_ 11« \ 

180 J 


76938698802699r 

12527254840352 


6211549394313 vr 3 
860382166336 


20491256151r 3 

7918618736 


186196891?; r 


2 2 


1778984341 v 3 r 77947-Cr __ 65875 vr 

2210193216 

( 12528089829r 4 
v 1124897312 


75 r 
1582 


40043584 
454075-id 2032591?; 3 


1236331119ur 

346090976 


lt _ 4 + 

1463625r 2 

2502724 


U- 3 + 


44296 69608 

7466191437-ur 3 22282569r 3 


646417856 


3652264 


111033673-ti 3 r 

66975552 


49 v r 
66 


15479 vr 
6328 


75r 

308 


4148928 


484 


2401i; 

39366 


24004512 

18974931^r 2 


2860256 


267775f 3 

889056 


49tT 

162 


1025-u 
3528 

, 103671r 2 
“ r 47432 

_ 25iA 


25v \ 
252 


21 1569041 160693r 4 _ 
12527254840352 

( | 3480659477t; 3 r 
l ' 2210193216 


7860884553897vr 3 , 85890653865r 3 , 186196891v 2 r 2 


84 J 

2047213695ur : 


860382166336 1 7918618736 

77947-u 2 r 198763^r , 2637r 


1 40043584 1 346090976 

454075-u 4 9562127u 3 . 1025^ 2 


44296 


69608 


1582 


4148928 


24004512 


3528 


J U-3 + 
U-2 + 

U_ 2 + 

13246263r 2 
2502724 

U-i + 


U - 1 + 


221v 

252 


28605737830234834755r' 

1467292304940749056 


2780975270379 15r 3 
21079363075232 


72837778155t> r 


2 2 


9690547328 

19447891v 2 r 50889r , 2706289217v 4 55091u 2 


8619732963r 

1211318416 


21 1569041 160693r 
12527254840352 


^ 1 5847072 

7860884553897i>r 3 


17402 


12098274048 


85890653865r 


186196891?; r 


63504 
2 2 


+ 1 ) liQ + 


2047213695^ 


U 0 


13246263 r 2 


860382166336 


7918618736 


40043584 


346090976 


3480659477t> 3 r _ 
2210193216 

( 12528089829r 4 


77947?; r 
44296 


198763t?r 

69608 


2637r 

1582 


454075?; 

4148928 


7466191437'ur 3 22282569r 3 


1124897312 
( i 111033673t> 3 r 


646417856 


3652264 


66975552 


49 v r 
66 


15479?;r 

6328 


75 r 
308 


484 

2401d 4 

39366 


9562127C 

24004512 

18974931i’r 2 

2860256 


1025-ti 

3528 


221v 

252 


2502724 
Mi - 


Ml 


1036711-- 

47432 


M 2 


_ 267775-» 3 , 49v 2 , 25t>\ , 

889056 ' 162 84 J “ 2 


76938698802699r 

12527254840352 


62 1 15493943 13iir 3 
860382166336 


20491256151r 3 

7918618736 


186196891?; r 


40043584 


162 


12363311 19 cd 
346090976 


1463625r 

2502724 


U3 


1778984341?> 3 r , 77947u 2 r , 65875-ur 


^ 2210193216 1 44296 1 

( 6296179542799261293r 4 _ 98816474558753811?>r 3 
2934584609881498112 


69608 


75r 

1582 


454075?; 4 


2032591?; 3 1025-u 2 


4148928 


24004512 


25-u 

252 


18549839506204160 


31068146188125r J 

42158726150464 


55178557611?; r 


3528 
287760797433tir 


19381094656 


133245025760 


401096399245-idr 

583491009024 


3588625 ur 
3898048 


573327 vr 
1093840 


75r 

8701 


1230500161?; 

24196548096 


268897» 3 

1714608 


16675?; 2 

127008 


M 3 + 

, 324338625r 2 

1 " 2422636832 

m) «•« 


The linearized viscous Burger’s equation (Eq. 15) is solved with this Compact scheme with multiple 
domain sizes to demonstrate the dependence of the implicitly derived spatial derivatives (and the stability 
limit) on the size of the domain. A Fourier stability analysis is performed and the stability of the Compact 
scheme as a function of Courant (r = ^r), Von Neumann (v = numbers, and (3 = Axk is derived as: 


R(G) = 


286057378302348347551- 4 _ 278097527037915r 3 __ 72837778155?; 2 ?- 2 , 8619732963r 2 

1467292304940749056 21079363075232 9690547328 ” l ~ 1211318416 

i 19447891ti 2 r _ 50889r , 2706289217?; 4 _ 55091?; 2 , 

5847072 17402 ~ l " 12098274048 63504 ~ l " 


(- 


21 1569041 160693r 4 , 85890653865r 3 , 186196891« 2 r 2 

6263627420176 _t ’ 3959309368 ' 20021792 


13246263r 2 \ 
1251362 J 


cos(P) + 


(- 


77947-» 2 r 

22148 


2637r 

791 


454075-u 4 

2074464 


1025t> 2 \ 
1764 J 


cos(/3) + 


( 12528089829)- 4 
^ 562448656 


22282569r 3 

1826132 


441-u 2 r 2 

242 


103671r 2 \ 
23716 J 


cos(2/3) + 


(- 


49 v 2 r 
33 


75 r 
154 


2401t> 4 

19683 


^ir) cos(2/3) + 
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M4 + 

(35) 



186196891t? 2 r 2 


1463625r' 

1251362 


(- 


76938698802699r 4 , 20491256151r 3 

6263627420176 ' 3959309368 20021792 


(+ 


77947 v 2 r 
22148 


75 r 

791 


454075m 4 

2074464 


1025m 2 \ 
1764 ) 


COS (3/3) + 


cos(3/3) + 


( 6296179542799261293)- 4 _ 31068146188125m 3 , 5517855761U’ 2 r 2 , 324338625r 2 \ 

1467292304940749056 21079363075232 ' 9690547328 ' 1211318416 J 


(- 


3588625m 2 m 

1949024 


150r 

8701 


1230500161m 4 

12098274048 


16675m 2 A 
63504 J 


COS (4/3) + 1 


COS (4/3) 


( 36 ) 


The imaginary part of the amplification factor is: 


3(G) 


+ 


( 


7860884553897ur 3 _ 2047213695 mm 2 _ 3480659477m 3 m 
430191083168 173045488 1105096608 



7466191437mm 3 

323208928 


18974931mm 2 

1430128 


111033673m 3 m 

33487776 


198763mm , 9562127m 3 
34804 12002256 


221m \ 

126 J 


sin(/3) 


15479mm _ 267775m 3 
3164 444528 



sin(2/3) 


( 62 115493943 13 ur 3 
^ 430191083168 


1236331119mm 2 

173045488 


1778984341m 3 m , 65875mm , 2032591m 3 

1105096608 + 34804 4 12002256 


25m\ 

126 J 


sin(3/3) 


98816474558753811mm 3 

9274919753102080 


287760797433mm 2 _ 401096399245m 3 m _ 573327 vr , 268897m 3 
66622512880 291745504512 546920 857304 



(37) 


The amplification factor is simply, |G| = y/3 ?(G) 2 + 3(G) 2 , and the range of stability (green area) is 
shown in Fig. 3. Notice how the stable region changes as the domain size, maxi (the number of grid points in 
each direction), changes. Explicit spatial derivative operators do not exhibit this behavior since the stencil 
size remains constant regardless of domain size. 


V. Application of UHF Method 

The MESA and Ultra Hi-Fi Methods are actually a procedure for designing ever more accurate numerical 
methods in which additional information is stored at each cell or grid point. For this comparison, the solution 
variable and up to it’s third spatial derivative will be stored at each grid. The notation, c4od, represents an 
UHF method with a 4 point stencil and only the solution variable and up to d derivatives on the grid. The 
basic procedure has been previously published for inviscid problems. 20 

The c4o0 UHF method will use a 4 point interpolation stencil to determine spatial derivatives as shown 
in Fig 4. A simple 1 st order Taylor series in time is used: 

< +1 = u n + u t (At) (38) 

in which the time derivative is found from the governing equation (Eq. 15) to be: 

du . . 

— = - cu x + im xx (39) 

and the solution variables, u, u x , u xx , are interpolated to the center of the four-point stencil using: 


< = Jq ("A-! + 9u i-l + 9u i+± ~~ u i + 1 
—u, a + 27u,_i ) — 27 u , , i + u t , s 

Tl L 2 2 ' 2 ' 2 

Xi ~ 24Aa; 








2Ax 2 


This results in a single time equation: 


(40) 

(41) 

(42) 


, n+U/2) 

a i 


/jAt cAt 1 \ 

~2h? 2Ah 16 y UUi ~i 


(iAt 9cAt 9 \ 

~ 2 h 2 + ~ 8 ^ + 16 j UUi ~i 


fiA t 9cAt 9 \ ( /iAt cAt 1 \ 

W ~8h~ + 16 ) uu ^ + V 2^ + m - 16 ) UUi +i 


(43) 
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(a) MAXI=2 


(b) MAXI=3 



(c) MAXI=4 


Figure 3. Lele Compact Linear Viscous Burgers Equation Stability Range 
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Figure 4. UHF Staggered Grid Diagram 


And after two time steps, the original ’’unstaggered” grid has been updated with the following: 
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(44) 


As was done with the Compact algorithm, a Fourier stability analysis is completed with the real part of 
amplification factor: 

3?(G) = (1152)- 1 (-2920u 2 + 288r(4r - 5) + 9 (348w 2 + 16(1 - 4r)r + 63) cos (/?) 

-18 (I2u 2 + 16r(4r - 5) + 9) cos(2/3) + (9(1 - 8r) 2 + 4c 2 ) cos(3/3) + 738) (45) 

And the imaginary part of the amplification factor: 

3(G) = t^ v (~ 8r + (8 r — 1) cos(/3) + 5)(sin(2/3) — 26sin(/3)) (46) 

48 

The amplification factor, |G| = \/3?(G) 2 + 3(G) 2 , is plotted as in Fig. 5 as a function of Courant and 
Von Neumann numbers, and (3 = Axk. The left figure (a) unwraps the polar plot in right figure (b). As 
(3 increases the amplification factor simply repeats. It has been traditional to use a polar plot due to this 
property. However, as soon as derivative information is stored on the grid, the amplification factor continues 
to vary as (3 > 2i r. This is due to the ability of these techniques to carry ultra-short wave information. 
Amplification factors larger than one for any wavenumber imply it is an unstable method. 

Following this same procedure, but adding additional derivatives on the grid to produce methods, c4ol, 
c4o2, and c4o3 results in the algorithms shown in the appendix. 
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(a) Amplification - c4o0 


(b) Polar Diagram 


Figure 5. c4o0 - Linear Viscous Burgers Equation 


The full equations shown in the appendix are required to perform the Fourier stability analysis to enable 
true method comparisons. The exact amplification factor is also shown since we know the exact solution. 
Notice in Fig. 6 that the various UHF methods more closely approach the exact ampification factor as the 
number of solution derivatives on the grid increases. Since the UHF methods can resolve ultra-short waves, 
the wavenumber range in the figure could be extended past k = 2ir. 

The stability region for all the UHF methods is shown in Fig. 7. This gives an indication of the allowable 
time step as the dynamic viscosity and convective terms vary. Notice (see definition of Von Neuman number 
given earlier) the role viscosity plays in reducing time step size. Fortunately, the viscosity, fi is a small term 
generally compared to convection, c. The Compact scheme has larger allowable time steps as shown in Fig. 3 
for a given grid. However, more grid is required for the Compact scheme as shown in the next section and 
despite this apparent advantage, a coarser grid actually results in an effectively larger time step for the UHF 
methods. 


VI. Application of the SIMPLE/PISO Methods 

The Semi-Implicit Method for Pressure-Linked Equations (SIMPLE) method 28 and Pressure Implicit 
with Splitting of Operators (PISO) method 29,30 are the mainstay of commercial fluid dynamics solvers. 

An overview of the SIMPLE method is: 

• Start the iterative process by guessing the pressure field. 

• Use those pressure values to determine the velocity from the momentum equations. 

• Determine a pressure correction such that the continuity equation is satisfied. 

• Find corresponding velocity corrections, and use new pressure and velocity. 

• Repeat this until a velocity field is found that does satisfy continuity. 

The governing equations are linearized to produce a system of linear equations with one equation for 
each cell in the domain. A point implicit (Gauss-Seidel) linear equation solver is used in conjunction with 
an algebraic multigrid method to solve the resultant scalar system of equations. 

The time stepping is first order accurate implicit, and the spatial accuracy is second order. 
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(a) r=v=.01 


(b) r=v=.05 



(c) r=v=.10 


(d) r=v=.5 


Figure 6. Ampification/Stability Factor Comparison 
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(a) Segregated Process 


Figure 8. Overview of Segregated Solution Technique 


The PISO algorithm moves the repeated calculations required by SIMPLE inside the solution stage of 
the pressure-correction equation to more closely satisfy the continuity and momentum equations. The PISO 
method takes more time per iteration, but often requires fewer iterations, particularly for transient problems 
as will be demonstrated next. 

Finally, the stability limits of these two approaches are large (limited by the need for accuracy) compared 
to the Compact and UHF methods due to implicit time-stepping. 
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VII. Nonlinear Navier Stokes — ID 


We will now test these techniques on an example of the one-dimensional Navier-Stokes equations that 
reduces to a one-dimensional heat transfer problem governed by the linear viscous Burger’s equation. This 
is what is solved in Sage, HFAST, and GLIMPS to produce reasonable Stirling analysis solutions. 

The one-dimensional Navier-Stokes equations can be written as: 23 


dU dE _ dE v 
dt + dx dx 


(47) 



pu 

pu 2 + p 

(E t T p)u 


E v 


0 

7~XX 

'U-Txx Qx 


E t = p 



T~xx 


4 

3 MW* 




k 


pCp 

Pr 


Sage/Glimps 32 utilize a slightly different form: 

dU dE _ dE v 
dt + dx dx 


(48) 

(49) 


(50) 



pu 

pu 2 + p 
(■ E t +p)u 


E„ 


0 

0 

Qi 


E t = p 




7~xx 


4 


<7; = - <k> 


dT p\u\ 
dx < p\u\ > 



(51) 

(52) 


Notice that Stoke’s stress tensor is replaced by source terms. This also simplifies the numerics since only 
a single second derivative, must be calculated compared to three in Eq. 47. We will utilize the standard 
Navier-Stokes form for comparison with commercial software. Extension of these ideas to the GLIMPS form 
is direct. 


A. Convection and Diffusion 

Few exact solutions exist for the Navier-Stokes equations and therefore validating commercial codes is typi- 
cally done with only approximate experimental information. One special case that both has an exact solution 
and yet includes heat transfer physics relevant to oscillating Stirling engines will be shown. 

First, commercial codes operate either in two or three dimensions. However, it is possible to cajole the 
commercial code into solving a one-dimensional problem by solving the full Navier-Stokes equations with 
an inherently one-dimensional problem such as flow through a pipe with an initial temperature ’’shock” as 
shown in Fig. 9. 

600K 300K 

u=l m/s — | — »- | u=l 

-2 m 0 m 2 m 

Figure 9. Heat Transfer Test 


This problem can be solved with either two or three-dimensional solvers and it reduces to the one- 
dimensional Navier-Stokes Eq. 47. Since the density and velocity are constant, both the continuity and 
momentum equations are satisfied. Only the energy equation actually needs to be solved: 


m 

dt + 


d_ 

dx 


^( pC v T + p)u 


-puu x + q x = 0 


(53) 
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Using, E t = pC v T , and dropping out terms that are zero, we have the following reduced form of the 
energy equation: 


_ W „ ST ,<PT 

= ‘a? (54) 

Finally, dividing by pC v and for essentially incompressible flows replacing C v with C p , we have the 
following transport (linear viscous Burger’s equation): 


dT dT d 2 T 
~dt +U lhc ~ “fV 2 


(55) 


with the thermal diffusivity a = . 

Furthermore, commercial solvers utilize the dimensional form of the Navier-Stokes, but by utlizing the 
following non-dimensionalization (* = nondimensional quantity): 


T* = — 

T 0 ’ 


•T ,* ^ * To * To 

= 7 , t = — ,u = u—, a = —ra 
L t 0 L L~ 


and using the following nondimensional boundary conditions: 


(56) 


T*(a:*,0) = 2 or 1, x € [-2,2], t* = [0,l], u* = 1.0, Ax* = .1, At* = .01, a* = .03003 
We have an exact solution, by the separation-of- variables, 33 (as N is large): 


N 


(2 k - 1) 


ir(x* - u*t*) 


L 


T*(x*,t*) = 1 + 0.5- - V sin 
7 r 

fc= l 

And finally, by choosing the characteristic constants as: 

To = Is, L = 1 m, To = 300/\ 


exp[— a*(2k — l) 2 7r 2 f*/T 2 ] 
2k- 1 


(57) 


(58) 


(59) 


We have the following problem definition for the commercial solvers: 

x £ [—2m, 2m], T £ [300AT, 600A'], t £ [0s, Is], u = 1 m/s, Ax = .lm, At = .01s, 

a = .03003m 2 /s, p = 998.2kg/m 2 , C p = 4182J/(kg • K), k = 125359J/(s • K), p = 101325Pa (60) 

This problem is solved and the results are compared as various techniques are applied to the problem. 
In Fig. 10, the exact solution is shown along with four UHF techniques, a 6 th order Compact scheme, a 
segregated spatially and temporally implicit method (SIMPLE) used in Fluent, a coupled spatially and 
temporally implicit method (PISO) used in Fluent, and a segregated central difference with a blended 
(averaged) Euler and Crank-Nicolson technique used in CFD-ACE. Fluent only allows 1 st order accuracy in 
time when moving meshes are applied as when the Stirling engine is simulated. However, CFD-ACE allows 
for up to 2 nd order accuracy in time when the meshes are compressed with a spring analogy. 

The best techniques in this test case are the Compact and c4o3 UHF schemes. The commercially used 
solvers are noticeably less able to model time accurate heat transfer. However, the coupled solver should 
be used when commercial codes are applied to oscillating Stirling simulations and when available, 2 nd order 
time accuracy should be used (as when using CFD-ACE). 

VIII. Fidelity and Turbulence Transition 

Modeling turbulence transition is a difficult problem due to the large disparity in both spatial and 
temporal scales caused when velocity gradients are high. In the Stirling engine velocity gradients are high 
near walls and regions of sheared flow due to oscillating/reversing flows. As the velocity gradients increase, 
the flow becomes rotational, leading to a vigorous stretching of vortex lines, which cannot be supported in 
two dimensions. 24 For this reason, truly turbulent simulations cannot be done in one-dimension. 
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Figure 10. Heat Diffusion and Convection Comparison 
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The Stirling engine exhibits turbulent and laminar behaviour simultaneously. 34 It is desirable therefore 
to avoid a Reynolds Averaged approach since the time averaged equations combined with some turbulence 
model 35, 36, 3 ' assumes turbulence everywhere. One promising approach is Large Eddy Simulation (LES) in 
which the Navier-Stokes equations are solved in time, but with spatial filtering applied, leaving the small 
eddies still unresolved. Since small eddies are essentially isotropic, the modeling is much easier compared 
to Reynolds time averaging. Moreover, the entire flow is bounded by walls making boundary condition 
specification much easier than the typical open domain problems encountered in modeling, for example, jet 
flow turbulence. 

The smallest scales of turbulence are the Kolmogorov scales of length, time and velocity: 38 

v = (v 3 / e ) 1/4 , t = (V/e) 1//2 , v = (i/e) 1 / 4 (61) 

where v is the kinematic viscosity, and e is the dissipation rate. The Reynolds number in the Kolmogorov 
region, Re = vr)/v = 1, shows the ratio of inertial and visous forces is unity because most of the energy is 
dissipated in this wavenumber region. 

In small Stirling engines we can estime the smallest scales as follows. With the average flow assumptions 
for Helium in the engine: 

g = 350e — 7(Ns/m 2 ), v = 502e - 6(m 2 /s), k = 278e — 3W/(mK), a = 29.9341e - 6m 2 /s, 

Pr = 0.654, oj = 2tt/ = 502.655(rad/s), C p = 5.19e3J/(kgK),T = 700/v, p = 2.6e6Pa, p = 1.78888(k g/m 3 ) (62) 

Then from West, 39 for oscillating flow, the average thermal boundary layer thickness, \J 2a jio = .345115mm, 
and the average flow boundary layer thickness, yj2vjui = 1.41329mm. 

The maximum surface shear stress may be approximated with this information by (assuming average 
maximum flow speed of 10 m/s and using the flow boundary layer thickness): 

A |/ 

T w = — = .247649A7m 2 (63) 

Ay 

The friction velocity, u T = = .372125m/s, is then used in the following equation for dissipation in 

channel flow: 38 

e « 2u 2 T U m /. 0253007m) = 54.7325m 2 /s 3 (64) 

Finally, we can estimate the Kolmogorv spatial wave length as: 

77 = (v 3 /e) 1 / 4 = 1.23301mm (65) 

Assuming a representative engine length of (1 in. = 25.3007?7im), the domain consists of 20 to 30 Kol- 
mogorov wavelengths, or roughly 8000 regions (Kolmogorov boxes) in 3D where isotropic turbulence mod- 
eling can be employed. This corresponds to a turbulent Reynolds number of Re T = u T (. 0253007m)/ (2^) = 

9.37751. This is low Reynolds number flow and appears ideal for applying Large Eddy Simulation. 

One would like to model the smallest turbulent scales, which in the case of LES is the Kolmogorov 
wavenumber range. The allowable time should also be of the same order as the Kolmogorov time scale. 

r = (v/e) 1 ' 2 = .00302s (66) 

This time scale is not prohibitive, since the time step size is already smaller than this in current commercial 
simulations due to the numerical issues involving moving grids. 

We would like to utilize the most efficient technique to minimize the computational cost. We know 
the minimum wavelengths required to simulate turbulence in the Stirling engine and the wave convec- 
tion/dissipation problem solution from Eq. 17 is used to determine the fewest grid points required for each 
method. The results are shown in Fig. 11. 

This solution represents a traveling (convecting) wave with a dissipating amplitude. Many methods will 
convect at the wrong speed (dispersive error) or will excessively dissipate the amplitude (dissipative error). 

A key measure of efficiency is how many grid points are required per wavelength to propagate this wave 
within some predetermined error bound. As k, the wavenumber in Eq. 17 , is increased, the number of grid 
points must increase to properly simulate the wave. 
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(a) Fidelity Comparison 


Method 

Spacing 

Error 

c4o0 

.1 

2.54229 10~ 2 

c4o0 

.2 

4.27563 10~ 2 

c4o0 

.4 

4.68577 KT 2 

c4ol 

.1 

3.11163 KT 6 

c4ol 

.2 

2.96551 10“ 5 

c4ol 

.4 

8.4702 10" 4 

c4o2 

.1 

1.0178 10" 10 

c4o2 

.2 

3.1935 10 -9 

c4o2 

.4 

6.41079 10" 8 

c4o3 

.1 

3.44169 10" 15 

c4o3 

.2 

2.27818 10" 13 

c4o3 

.4 

3.12925 10" 11 

Compact 

.1 

1.0993 10 -6 

Compact 

.2 

7.10929 10" 5 

Compact 

.4 

5.31245 KT 3 


(b) Errors at t=l 


Figure 11. Wave Propagation Fidelity 
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Currently, Compact schemes are regarded as requiring approximately 6 grid points per wavelength for 
reasonable solutions. A comparison of the Compact scheme and the UHF schemes with up to 3 spatial 
derivatives stored on the grid, as shown in the table (Fig. 11(b), demonstrates the Compact scheme is 
comparable to method c4ol (one solution derivative per grid point). This is expected since the Compact 
scheme also utilizes 1st derivative information. However, extending Compact schemes to include higher 
derivatives involves complicated matrices which may be difficult or intractable to solve. 

The c4o3 method can match the results of the Compact scheme using 16 times fewer grid points per 
dimension. Specifically, one additional test of the c4o3 method with k = 8 tt and Ax = .2, the error at t, = 1 
was 2.44291 5, which is more accurate than the Compact scheme with considerably fewer grid points per 

wavelength. Note that the previous comparisons had k = tt and here it was multiplied by eight and we use 
half as many grid points for a net change of sixteen grid points per wavelength. 

Clearly, the UHF schemes are more efficient and can be formulated explicitly for easier parallelization. 
This implies smoothly transitioning turbulent flows can be more efficiently simulated with UHF techniques. 

IX. Conclusion 

One-dimensional Navier-Stokes equations are currently utilized for Stirling engine design and optimization 
with reasonable success. Recent attempts at multi-dimensional simulations have relied upon commercial 
solvers and this report examined that practice more closely. 

This report has shown that the techniques used in commercial codes for simulating Stirling engines are 
not as capable as more recently developed approaches available in the literature. Moreover, the unique 
environment of the Stirling engine in which flow is simultaneously turbulent and laminar makes large eddy 
simulation desirable, while the low Reynold’s number, wall bounded flow provides for modest grid require- 
ments and well defined boundary conditions. 

Despite the larger stability limit (4 times larger) of Compact schemes for a given grid spacing, the UHF 
method results in an effective time-step that is 4 times larger than Compact schemes since the grid can be 16 
times coarser per dimension. The 6 th order Compact schemes performed well with the heat transfer test and 
apparently would work well in regions of conjugate heat transfer. However, the Compact scheme is not as 
efficient at predicting turbulent transition compared to UHF methods. The c4o3 method performs similarly 
to the Compact scheme for steep temperature gradients (conjugate heat transfer) but is up to 16 3 = 4096 
times more efficient when three-dimensional transitional flows need modeled. 

It would be desirable to compare c4o4 methods and higher in the future. Future work should exam- 
ine utilizing UHF methods in a steady-harmonic formulation with Detached Eddy Simulation and more 
complicated moving grid tests should be performed. 
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X. Appendix 


In what follows are the full equations used to time advance the solutions in this paper and their stability 
amplification factors. These equations are lengthy, but are provided for completeness and to allow for 
independent verification. 

A. c4ol 

The single step c4ol is given by: 
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256 


/ 70m 3 At 3 
^ 3ft® 


15cm 2 Ai 3 55c 2 mA* 3 , llc 3 At 3 55 m 2 A t 2 , IIcmA* 2 , 19c 2 Ai 2 , 19 mA* 

ft 5 12ft 4 + 24ft 3 12ft 4 4ft 3 + 192ft 2 ^ 96ft 2 


19cAi 

192ft 



( 68 ) 


The c4ol full step (two staggered half-steps) is: 


n+1 


(69) 


/ 770m 3 At 3 | 45 cm 2 Ai 3 785 c 2 mA* 3 15c 3 Ai 3 785 m 2 At 2 45 cmA* 2 , 281c 2 At 2 , 281 mA* , 27cAt 5 \ 

9ft' 3 ft® 36ft® 8ft 4 36ft 5 ift 4 576h. s ' + 288ft 3 3 " 64ft 2 256ft ) “*-§ 


+ 


( 


210m 3 Ai 3 
ft 7 


45cm 2 Ai 3 

P 


345c 2 mA* 3 

4ft 5 


15c 3 A* 3 345m 2 Ai 2 . 45cmA* 2 

8ft 4 4ft 5 3 ~ 4ft 4 


513c 2 At 2 , 513 mA* 27cAi 405 \ 

64ft :! ' 32ft 3 64ft 2 256ft y “i-j 


+ 

+ 


(- 

(- 


210m 3 Ai 3 45cm 2 At 3 

ft 7 ft® 

+ 

345c 2 mA* 3 , 

4ft 5 3 " 

15c 3 Ai 3 
8 ft 4 

+ 

345/i 2 A/ 2 | 45c/iA/ 2 
4ft 5 1 4ft 4 

513c 2 At 2 
64ft 3 

513/iAt 

32ft 3 

27cAt 

64ft 2 

+ 

405 ^ 
256ft y 

K+i 

770m 3 At 3 i 45 cm 2 Ai 3 

+ 

785c 2 mA i 3 

15c 3 A* 3 

+ 

785m 2 At 2 

45c/xAt 2 

281c 2 At 2 

281mA* , 

27cAt 

+ 

5 > 

+ i+ | 

9ft 7 1 ft® 

36ft 5 

8 ft 4 

36ft 5 

4ft 4 

576 ft 3 

288ft 3 3 ~ 

64ft 2 

256ft J 
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, ( 70/i 3 At 3 
+ ( 3 ft 6 

, ( 210/t 3 A t 3 

+ ( ft 6 

, ( 210/i 3 At 3 

+ I ft 6 


l 3/i 6 


, 15c/i 2 At 3 

+ ft 5 

45c/i 2 At 3 

! w 

45c/? 2 At 3 
ft 5 


ft 5 


55c 2 ^At 3 


11c 3 Ai 3 

55/t 2 A t 2 

llc/?Ai 2 , 

19c 2 A? 2 , 19/iAt , 

19cAt 

256 ) 

12 ft 4 


24 ft 3 

12 h 4 

4 h 3 

192ft 2 1 96ft 2 ” l " 

192 h 

285c 2 /?At 3 


19c 3 At 3 

285/i 2 At 2 

57cfiAt 2 , 

297c 2 At 2 , 297a?A t 

1 99cA£ 

_ 81 ^ 

4 ft 4 


8ft 3 

4 h 4 

4 h 3 

64ft 2 1 32ft 2 

64 h 

256 y 

285c 2 /?At 3 

+ 

19c 3 At 3 

285/? 2 A t 2 , 

57c/j,At 2 ■ 

297c 2 At 2 , 297a?A t 

99cAt 

_ S' 1 ) 

4 h* 

8 h 3 

4ft 4 

4/? 3 _l_ 

64ft 2 1 32ft 2 

64 h 

256 J 

55c 2 fiAt 3 

+ 

11c 3 At 3 

55/? 2 At 2 , 

llc/?At 2 , 

19c 2 At 2 , 19/xAt 

19cAt 

256 ) U - 

12 ft 4 

24ft 3 

12ft 4 

4ft 3 + 

192ft 2 1 96ft 2 

192 h 


i- a 


u x 


i+i 


(70) 


The real part of c4ol amplification factor: 


sft(G) = (13824 -1 )(32 (1084?; 3 — 3(8r(1160r — 271) + 225)i> + 216) cos(/3) + 

(2248n 3 + 972(80r - 3)v 2 - 48r(6280r - 281)?; - 270?; - 648r(40r(8r - 3) + 9) + 35l) cos(2/3) - 
48/3 (5760r 3 + 48(400?; - 51)r 2 - 6 (408t> 2 + 436?; -79 )r + n((237 - 436?;)?; + 90)) sin(/3) + 

81 (2560r 3 + 320(46?; - 3)r 2 - 24(2?>(20?; + 57) - 3)r + Qv ( -76w 2 + 6u + 45) + 28/3sin(/3) + 8l) - 
3/3 (23040r 3 + 2112(10?; - 3)r 2 - 24 (264n 2 + 38?; - 19) r - (2v - 3) (76?; 2 - 9)) sin(2/3)) (71) 

The imaginary part of c4ol amplification factor: 


3(G) = (13824 _1 )(— 81/3 (2560r 3 - 1216(10?; + l)r 2 - 8(2?/(76?; - 99) - 33)r + 3(2?; + 1) (44?; 2 - 9)) + 

6/3 (9(82?; - 39) + 4 (-910?; 3 - 3(1104r - 139)?; 2 + 5460r(8r - 1)?; + 6r(24r(40r - 23) + 139))) cos(/3) + 
3/3 (23040r 3 + 2112(10?; - 3)r 2 - 24 (264?; 2 + 38?; - 19) r - {2v - 3) (76?; 2 - 9)) cos(2/3) 

+2 (19592?; 3 - 972(80r - 3)?> 2 - 6(8r(15560r - 2449) + 1845)?; + 27(24r(40r(8r - 3) + 9) + 115)) sin(/3) 
+ (2248?> 3 + 972(80r - 3)?; 2 - 48r(6280r - 281)?; - 270?; - 648r(40r(8r - 3) + 9) + 35l) sin (2/3)) (72) 

B. c4o2 

The single step c4o2 (with r = n-^s and v = | is: 


,"+( 1 / 2 ) 

a i 


(oooc.,5 29575«r 4 7875 r 4 7875u 2 r 3 , 113995 vr 3 \ , 

4 8 4 1 144 ) U i~i + 


4425 r 3 , 113995v 3 r 2 , 13275?; 2 ?' 


( 


32 1 

I 4425 v 4 r 
' 128 

22223t> 5 _ 
41472 


288 1 
111115-u 3 r 


64 

2025u 2 r 


111115i>r 

3456 

13297vr 


_ 2025r 
256 

1215r\ 


(- 
(- 

( ooocJ I 250425vr 4 , 7875r 4 , 7875 v 2 r 3 148365-ur 3 4425r 3 A , 

H 2 I g I 2 16 32 J u i- 5 + 


10368 256 1 18432 1 4096 j 

675?; 4 , 13297? 3 , 1215?; 2 155?; 383 


1024 


110592 


8192 


16384 


32768 


(' 


(- 


148365v c r 


13275v r 


98415ur 


32 64 1 128 

2025r 2 4425v 4 r , 32805v 3 r , 2025v 2 r 80919ur 1215r 


256 

I 6561-c 5 
r 512 


128 

675t> 4 

1024 


128 

26973t? 3 

4096 


256 

1215t; 2 

8192 


2048 

32805f j 
16384 ' 

,,2 3 


Uj_ i + 


4096 

16767 A , 

32768 J ^ 


— 2835r 5 — 250425?;r 4 _j_ 7875 r 4 _|_ 7875v 2 r 3 _|_ 148365?;? 
* 4 8 4 


16 


Ua 


4425 r 3 
32 


148365-iGr 

32 


13275v 2 r 2 

64 


98415fr 2 

128 


?+i + 
'*+4 + 


( , 2025t£ _ 4425?; 4 r _ 32805?; 3 ?- , 2025?; 2 r , 80919?;? _ 1215r\ 

256 128 128 256 2048 4096 J U ?+j 


(- 


6561?; 5 

512 


675u 4 

1024 


26973?> 3 

4096 


1215u 2 _ 32805?; , 16767 A 
8192 16384 32768 J “?+j 
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2835r 5 


29575t?r _ 7875 r 4 _ 7875iTr’* _ 113995^ , 4425r 3 

4 8 4 144 ' 32 


113995t; r 
288 


UaiS + 


13275t; 2 r 2 , 111115^r 2 _ 2025r 2 , 4425^ 4 r , 111115t; 3 r 

64 ' 3456 256 ' 128 ' 10368 


Uaa_ 3 + 


2025^ 2 r _ 13297fr , 1215r , 22223?; 5 _ 675t; 4 _ 13297^ 3 

256 18432 ' 4096 ' 41472 1024 110592 


110592 j a i+i ' 


16384 32768 ) u i+i 


Uaj ,3 + 


(l05/i 2 r 5 - - 7 f/, 2 rr' + - 2 -fh 2 v 2 r 3 + /i 2 w 3 ) u** + 

( | 805 1,2 3^2 149 h 2 r 2 . 345 ft 2-.2„2 745ft 2 ?.?- 2 \ , 

g6 II u ' 768 1- 64 n u ' 1152 y 4~ 

/ I 115 i'2 ? .4 r _ 745 h 2 v 3 r , 29ftV _ 149 t2 2„ , 29h 2 vr _ 149h 2 v 5 \ , 

\^ _r 128 1 u 1 3456 “ r 4096 768 n u ' " r 2048 13824 y a xx t _3 4~ 

/ 149ft 2 ?; 4 I 29ft 2 ?. 3 _ 9h 2 , 29ft 2 ?; 2 __ 3h 2 v \ , 

9216 4~ 12288 32768 4~ 8192 16384 J a xx i _3 ' 

(-945 h 2 r 5 + + 27^5^4 _ + ^Vr 3 ) + 

(_ 45 ^ 5^3 _ i^w + _ «™>fcVr 2 ) UxXi _ h + 

( 4 - ?315 z. 2 2 _ 2235/,2,,4 t , , 3105 1,2,, 3„ _ 1701ft 2 r , 1863;, 2 ,,2,. _ 5103ft 2 ?;?- , 621ft 2 ?. 5 \ , 

^4~ 128 1 128 n u ' 128 n u 4096 4~ 2 56 y ' 2048 4~ 512 J u xx i _i ~r 

( i 621 ft 2 ?; 4 _ 1701h 2 v 3 i 729 h 2 _ 1701ft 2 ?. 2 , 729ft 2 ?; \ , 

^4~ 1024 4096 4~ 32768 8192 4~ 16 384 y u xx i _i ~r 

(-945 /i 2 r 5 + _ 2^05 _ 2235|M + scms/^s + + u + 

/ I 1863ft 2 ?- 2 _ 6705 l2,,2„2 _ 9315 h 2,, r 2 __ 2235 ;,2„4 r _ 3105/,2„ 3„ _ 1701ft 2 r , 1863 l 2„2„ , 5103ft 2 ?;r- \ , 

^ 256 64 " u ' 128 11 ul 128 n u ' 128 u ' 4096 2 56 n u ' 2048 y Ua:a: i + i ^ 

/ 621/i 2 t> 5 I 621/i 2 u 4 I 1701/t 2 t> 3 , 729 h 2 _ 1701h 2 v 2 _ 729h 2 v \ , 

^ 512 + 1024 4096 327 68 8192 16384 J u xx i+ i ~r 

(l05 h 2 r 5 - 245|!rl + lf h 2 vr A + h 2 vr 3 - ff h 2 v 3 r 2 - + ^/i 2 n 2 r 2 +) u XXi+ + 

( 745h 2 vr 2 , 115 l 2.,4„ , 745/i 2 « 3 r , 29h 2 r _ 149 u2.,2 r _ 29 h 2 vr \ , 

^ 1152 4- 128 " t' ' 4- 3456 4- 4096 7 6 8 ,t ' 2048 y 4- 

/ I mV _ 149/i4V __ 29h 2 v 3 _ 9 h 2 , 29h 2 v 2 , 3h?v_\ , 

^ 13824 9216 12288 32768 8192 + 16 384 y U:Ex i + 3 4" 

(l015ftr 5 - - ^-hv 2 r 3 + + ^tor 3 ) w x ._3 + 

I 10535 l. 3 r 2 I 4055 u 2 2 _ 54011tr 2 _ 9995hvr 2 . \ , 

y 4” 96 ' 4- g 4 /tt' / 2304 1152 4-y 4- 

/ 4055 l 4„ _ 9995 hv 3 r _ 5401 hv 2 r , 355fer , 1181/tur _ 1999ft,» 5 _ 5401ft,t> 4 \ , 

^ 384 UL 1 3456 2304 4" 4096 4” 6144 13824 27648 J Ux i _ | 4" 

( | 1181/tf 3 355hv 2 _ lllh _ 41feu \ , 

1^4" 36864 4" 8192 32768 16384 y Ua; i-a 4" 

(945/ir 5 - 220 f r4 + §5775 hvr 4 _ m to 2 r 3 + 55||r^ _ 58^5^3^ ^ + 

(-»/n; 3 r 2 + iff to 2 r 2 + ^ tor 2 + f§^ 4 r) u Xi _ s + 

f 14715 u 3 , 1377; 2 _ 6075/ir _ 38637hvr , 2943/u; 5 , 459ftt> 4 __ 12879/t » 3 _ 6075hv 2 , 5103ft , 8019ft« A , 

^ 128 lU ' 256 4096 2048 “ r 512 1024 4096 8192 “ r 32768 _t ’ 16384 y 

( — 945/ir 5 + “QSfcrl + 89775 tor 4 + » to 2 r 3 _ 55|ftr^ _ 58695^3 _ 58695^2 _ 1665^2^^ + 


16384 J 3 


(-^pWV 2 + to 2 r 2 + W*._ i + 

1377 l 2 _ 6075/ir _ 38637 hvr , 2943/ii; 5 , 459h^ 4 _ 12879^ 3 _ 6075/i^ 2 , 5103/i , 8019 /it; 

256 ' 4096 2048 ' 512 ' 1024 4096 8192 ' 32768 ' 16384 


1377frr 2 i 44145 
256 ~ 1_ 128 


+ - m hv4r ) u x i+ i + 


3 _ 1377 u 2 6075ftr _ 38637hvr , 2943ft« 5 

256 4" 4096 2048 4" 512 


459hv 4 12879ft?; 3 , 6075 ft?; 2 5103ft , 8019ft?> 


32768 1 16384 y x i+i 


NASA/TM— 2005-213976 


23 



(-1015ft.r 5 


7945/ir 

24 


( 5401 hr 2 
\ 2304 


(- 


+ ^hv 2 r 3 - 405 ^ r3 + ^§^hvr 3 + hv 3 r 2 - 

9995hvr 2 4055 UnA r 9995 hv 3 r , 5401 hv 2 r 355 hr _i_ 1181 hvr 

1152 384 ,iU ' 3456 ' 2304 4096 ' 6144 

1999 hv 5 i 5401fof 4 , 1181hv 3 _ 355 hv 2 , lllh _ 41 hv \ 

13824 ' 27648 ' 36864 8192 ' 32768 16384 ) Ux i+% 


4055 l 2„2 
64 ,lU ' 


(73) 


The double step from it™ to w" +1 can be immediately derived by using two single steps above. The 
equations are ommitted due to length. 


1. c4o2 Amplification 

The equations for the amplification factor of this method are sufficiently large that their Fortran form is 
shown instead. 


\Im(G)= (beta*v* (-130029412509 - 135591927894717235200*r**8 + 

235688089173884928000*r**9 + 2053416202920*v**2 - 
12925760780304*v**4 + 37818812555520*v**6 - 
41659973309440*v**8 - 

59454259200*r**7* (-629819125 + 3450959988*v**2) + 
6193152000*r**6* (-1050833123 + 10984215328*v**2) + 

7372800*r**5* (105946873317 - 1571117362052*v**2 + 
2866080408576*v**4) - 

368640*r**4*( 179452656305 - 3380823300140*v**2 + 

1 1635255448992*v**4) - 

122880*r**3* (-33255849510 + 726770729579*v**2 - 
3540668800920*v**4 + 3497307329712*v**6) + 

192*r**2* (-955513044297 + 21845158240800*v**2 - 
132867547584800*v**4 + 234180801868800*v**6) + 

160*r* (37334893005 - 789506398980*v**2 + 5296053262104*v**4 - 
12801424523744*v**6 + 7884402819840*v**8) ) ) /7 . 3383542784el0\ 
+ (beta*v* (692362647 - 5462028784867737600*r**8 + 

946952243616 1536000*r* *9 + 54343258344*v**2 - 

53539971 1056*v**4 + 1707900809472*v**6 - 1864301675520*v**8 - 

138726604800*r**7* (-11074625 + 58913076*v**2) + 

2890137600*r**6* (-95039991 + 964377760*v**2) + 

491520*r**5* (69443249579 - 998375458220*v**2 + 
1771442856960*v**4) - 

40960*r**4* (72884695365 - 1333734152300*v**2 + 
4436424946464*v**4) - 

40960*r**3* (-4525470270 + 98467155981*v**2 - 
464360171848*v**4 + 442587839568*v**6) + 

64*r**2* (-118439206263 + 2973197651040*v**2 - 
17955426 117600*v**4 + 30755816238080*v**6) + 

32*r* (4848891741 - 163397666340*v**2 + 1199037538872*v**4 - 
288590 1281760*v**6 + 1738612309760*v**8) ) *Cos (beta) ) / 

3 . 623878656e9 - (beta*v* 

(8037225 - 747508888240 12800*r**8 + 129148920004608000*r**9 - 
34263000*v**2 - 575412336*v**4 + 5181034752*v**6 - 
105494108 16*v**8 - 69363302400*r**7* (-276055 + 1711536*v**2) + 
82575360*r**6* (-34655721 + 424199720*v**2) + 

245760*r**5* (1083499015 - 20285505016*v**2 + 

43052870784*v**4) - 

40960*r**4* (362625615 - 10128202520*v**2 + 45084965856*v**4) - 
2048*r**3* (-224343621 + 9745619820*v**2 - 72360806480*v**4 + 
92721195840*v**6) + 

64*r**2* (-76222053 + 7530750144*v**2 - 92593649088*v**4 + 
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236994878720*v**6) + 

8*r* (-22830093 - 572822928*v**2 + 13259861136*v**4 - 
59079642624*v**6 + 50087375360*v**8) ) *Cos (2*beta) ) / 
2.26492416e8 - (beta*v* 

(-164355237 + 534389864162918400*r**8 - 
2562002937446400000*r**9 + 3250233864*v**2 - 26360268048*v**4 + 
94274868480*v**6 - 127264420864*v**8 + 

59454259200*r**7* (607585 + 17586156*v**2) + 

247726080*r**6* (-81048203 + 261110560*v**2) - 

1843200*r**4* (130398037 - 2750451916*v**2 + 9383105312*v**4) + 

1474560*r**5* (1943667693 - 29362137364*v**2 + 

22507355136*v**4) - 

24576*r**3* (-538843482 + 12803987455*v**2 - 71553575640*v**4 + 
66154527600*v**6) + 

192*r**2* (-2579408793 + 65572339872*v**2 - 457570639904*v**4 + 
1041142379520*v**6) + 

32*r* (375564033 - 8984298708*v**2 + 70710058008*v**4 - 
228004837984* v**6 + 212775217920*v**8) ) *Cos (3*beta) ) / 

2 . 93534171 136e 11 + (v* 

(beta**2* (261842949 - 2194527814366003200*r**8 + 
3724642866954240000*r**9 + 6159355992*v**2 - 
953771 69520*v**4 + 415611350784*v**6 - 559729462272*v**8 - 
6606028800*r**7* (-91437875 + 508551372*v**2) + 
1513881600*r**6* (-67322291 + 726020640*v**2) - 
40960*r**4* (21099636639 - 445669712500*v**2 + 

16346 11927648*v**4) + 

163840*r**5* (70054445343 - 1103008857580*v**2 + 
2070815523840*v**4) - 

40960*r**3* (-1054627614 + 28532469783*v**2 - 
1564077 13592*v**4 + 164686577776*v**6) + 

64*r**2* (-19978382781 + 697506487200*v**2 - 
5261467218336*v**4 + 10338925880320*v**6) + 

32*r* (489426543 - 28430441 100*v**2 + 286355835144*v**4 - 
845957587872*v**6 + 574710288640*v**8) ) - 
16* (134555175 - 1115128841895936000*r**8 + 

1927070689001472000*r**9 + 9672938919*v**2 - 
77637693564*v**4 + 255245658912*v**6 - 304576780800*v**8 - 
111476736000*r**7* (-2734475 + 15275008*v**2) + 
1135411200*r**6* (-45562643 + 488318640*v**2) - 
307200*r**4* (1552020979 - 30919444147*v**2 + 
110645478832*v**4) + 

61440*r**5* (97159680473 - 1493878846040*v**2 + 
2789562424320*v**4) - 

20480*r**3* (-1325701170 + 31436422945*v**2 - 
162423477334*v**4 + 166958390760*v**6) + 

80*r**2* (-13555495233 + 349150858992*v**2 - 
2308931970240*v**4 + 4292918703616*v**6) + 

4*r* (6887181573 - 188364602520*v**2 + 1419968235096*v**4 - 
3707569366400*v**6 + 2391656357120*v**8)))*Sin(beta))/ 

3 . 623878656e9 - (v* (1741890648224563200*r**8 - 
3029855637602304000*r**9 + 

178362777600*r**7* (-2524059 + 15419236*v**2) - 
82575360*r**6* (-835546259 + 9961097280*v**2) + 

122880*r**4* (3233787985 - 85348974940*v**2 + 

3630679 15 168*v**4) - 
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49 1520*r**5* (13630408051 - 245209154932*v**2 + 
507285257472*v**4) + 

40960*r**3* (-324619695 + 13089759659*v**2 - 91073759432*v**4 + 
111515155248*v**6) - 

64*r**2* (-2872614537 + 218229409440*v**2 - 2469996553632*v**4 + 
59749 1 9813 120*v**6) + 

3* (-65098971 + 157579992*v**2 + 6234027408*v**4 - 
48443322624*v**6 + 92309828608*v**8) - 
64*r* (-53262927 - 2462239134*v**2 + 47336581836*v**4 - 
197637357392*v**6 + 159618062720*v**8) + 

2*beta**2* (4153113 - 37031679885312000*r**8 + 

73234329580339200*r**9 - 39990024*v**2 + 24132816*v**4 + 
1091715840*v**6 - 3011607552*v**8 - 
33030144000*r**7* (-249575 + 1831788*v**2) + 

27525 120*r**6* (-38982313 + 549851040*v**2) - 

40960*r**4* (91581147 - 3105008708*v**2 + 16241654240*v**4) + 

163840*r**5* (519280935 - 11052165404*v**2 + 

27241385472*v**4) - 

8192*r**3* (-7629228 + 607651455*v**2 - 5642488280*v**4 + 
8549724400*v**6) + 

64* r **2* (34845471 + 1103282208*v**2 - 23979806496*v**4 + 
71777561600*v**6) + 

32*r* (-5093523 + 25329348*v**2 + 600224904*v**4 - 
3713708832*v**6 + 3396972800*v**8) ) ) *Sin(2*beta) ) / 

1 . 811939328e9 - (v* (-3*beta**2* 

(4516155 - 67371 18835507200*r**8 + 57173118183014400*r**9 - 
90405720*v**2 + 745088112*v**4 - 2739588864*v**6 + 
3881749504*v**8 - 178362777600*r**7* (15625 + 90076*v**2) - 
247726080*r**6* (-3042061 + 20646880*v**2) - 
1474560*r**5* (62705127 - 1050002828*v**2 + 

1304454144*v**4) + 

368640*r**4* (19645223 - 434387796*v**2 + 1661380448*v**4) + 
24576*r**3* (-15672258 + 381876625*v**2 - 2275528200*v**4 + 
2407898640*v**6) - 

192*r**2* (-73413567 + 1893294048*v**2 - 13708750304*v**4 + 
32725847040*v**6) - 

32*r*(10491039 - 254143116*v**2 + 2039834376*v**4 - 
6769699744*v**6 + 6493098240*v**8) ) + 

4* (139405941 - 634572669085286400*r**8 + 

24298575227781 12000*r**9 - 2722769964*v**2 + 

21720929760*v**4 - 75428706432*v**6 + 96443866112*v**8 + 
1981808640*r**6* (5826243 + 14884490*v**2) - 
29727 129600*r**7* (-498163 + 39563832*v**2) - 
15482880*r**5* 

(131848619 - 1741393344*v**2 + 166073344*v**4) + 
2580480*r**4* (71794645 - 1435990910*v**2 + 

4239 126 144*v**4) + 

122880*r**3* (-86251770 + 1995280729*v**2 - 
10367448084*v**4 + 8023615488*v**6) - 
384*r**2* (-1055068551 + 26440604040*v**2 - 
1772903 19536*v**4 + 382393509120*v**6) - 
16*r* (6263857 17 - 14795770608*v**2 + 114091371864*v**4 - 
356757419584*v**6 + 323542759680*v**8)))*Sin(3*beta))/ 

2 . 93534171 136e 11 
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C. c4o3 

The single step c4o3 is : 


n+( 1/2) 


( 1 oci ocn„7 i 819043225vv 6 , 1819125v 6 , 5457375v 2 v 5 272189225vv 5 2163105v 5 1360946125v 3 v 4 A , 

^ l^oioour -I- 162 + 4 + 4 432 32 2592 ) U 'i~ | 


(- 

(- 


10815525v 2 r 4 

64 

1097084975vr 3 


432659675vr 

10368 

515025r 3 


1414875r 

256 


86531935v r 


3605175v r 
64 


1414875v r 


432659675v r 
15552 


1097084975v c V 


1414875v r 
128 


1545075v r 


746496 2048 1 20736 1 512 1492992 4096 

( i 7170265vr 2 , 127575r 2 , 94325v 6 r 219416995v 5 r 515025v 4 r , 7170265v 3 r 


( 


221184 1 16384 1 512 2985984 

i 127575v 2 r 286345vr 25515r 31345285v 7 


8192 1 663552 

34335v 6 , 1434053v 5 


Ua 


16384 589824 131072 17915904 16384 1 2654208 

42525v 4 286345v 3 25515v 2 , 7435v . 5969 


(- 


65536 


3538944 


262144 


1572864 


1048576 


U. 


i— A 


/ionocn„7 , 115540425vv 6 1819125v 6 5457375 v 2 v 5 144542475 vv 5 , 2163105v 5 240904125v 3 v 4 A , 

^1<K>130U7 i 2 4 4 16 1 32 32 J U i + J + 


(+ 

(- 


2 4 


10815525 vV 
64 

52367175vr 3 


108103275vr 


1414875r 


128 

515025r 3 


256 

21620655v 5 r 2 


3605175v r 
64 


1414875v r 


36034425v r 
64 


52367175v c V 


1414875v r 
128 


1545075v r 


H+ 


I + 


1024 1 2048 1 

16996635vr 2 127575r 2 


256 512 2048 1 4096 

94325v 6 r 10473435v 5 r , 515025v 4 r , 5665545v 3 r 


u i+k + 


( 


8192 16384 512 4096 

_127575v 2 r 4275585vr , 25515r 498735v 7 


1 8192 1 8192 

34335v 6 , 1133109v 5 


u i+b + 


16384 65536 1 131072 8192 1 16384 1 32768 

42525v 4 1425195v 3 , 25515v 2 , 1213785v , 518319 A 


h+i + 


(- 


65536 


131072 


262144 


524288 1 1048576 ) 


Ua_ i + 


( ioKiQCin.,7 115540425vr 6 1819125r 6 5457375v 2 r 5 , 144542475vr 5 , 2163105r 5 , 240904125v 3 r 4 A , 

^lOOlODUr 2 4 4 1 16 1 32 1 32 ) U i+\ + 


(+ 

(- 


2 4 


10815525v r' 
64 

, 52367175vr 3 


108103275vr 


1414875r 


128 

515025r 3 


256 

21620655v 5 r 2 


3605175v r 
64 


1414875v r 


36034425v r 
64 


52367175v’V 


1414875v r 
128 


1545075v r 


H+ 


I + 


1024 1 2048 

16996635vr 2 127575r 2 


256 512 1 2048 1 4096 

94325v 6 r , 10473435v 5 r , 515025v 4 r 5665545v 3 r 


u i+b + 


( 


8192 16384 512 1 4096 

127575v 2 r . 4275585vr . 25515r . 498735v 7 


1 8192 8192 

34335v 6 1133109v 5 


U i+b + 


16384 1 65536 1 131072 1 8192 1 16384 32768 

42525v 4 , 1425195v 3 , 25515v 2 1213785v , 518319 A 


h+i + 


( 


(- 


65536 


131072 


262144 


524288 


U iA-l + 


1048576 ) u i+? 

_1 or i QCfw-7 _ 819043225vr 6 , 1819125r 6 , 5457375v 2 r s , 272189225vr 5 _ 2163105r 5 , 1360946125v 3 r 4 \ , 

lOOlOOUl 162 ' 4 ' 4 ' 432 32 ' 2592 J 


(- 

( 4 


10815525v 2 v 4 

64 

1097084975vv 3 

746496 


432659675vv 4 

10368 

515025v 3 


1414875v 

256 


(- 


86531935v r 

2048 20736 

94325v 6 r 
512 ' 

25515r 

16384 "H 589824 131072 

42525v 4 i 286345v " 


3605175v r 
64 


1414875v r 
512 


432659675v r 
15552 


1097084975v ,j r 

1492992 


1414875v r 
128 


1545075v r 


4096 


U i+% + 


U i+ 3 + 


7170265 vr 2 , 127575r 2 
221184 ' 16384 


219416995v 5 r _ 515025v 4 r 
2985984 8192 


7170265v 3 r 

663552 




127575 v z r 


286345 vr 


(- 


31345285v 7 _ 34335v 6 _ 1434053v 5 
~~ r 17915904 16384 2654208 

25515v 2 7435v , 5969 


u i + 1 + 


65536 


3538944 


262144 


1572864 


1048576 


U. 


'i+¥ 


^_mmO h 2 r 7 + 63140h a v» + 1796795^2^6 _ 2313|ftV + QZWdhWr 5 - ^^h 2 Vr 5 - ^ f^h 2 V 3 r 4 ^) U XXi _ 3 + 

( | 188405/i 2 r 4 _ 115675 l 2„,2„4 , 818965 ^2,^4 _ 115675 t.2,,4^3 , 818965 L2,.3 y .3 _ 131695/i 2 r 3 \ , 

^ 864 16 11 U 1 " r 576 n Ul 48 U 1 " r 864 U U 1 13824 J a xx i _3 


(- 


188405 l 2 2 3 _ 2000635fr 2 vr 3 , 163793^ 2 v 5 r 2 
432 ,l U ' 41472 1152 


188405/i v r 
1728 


2000635/1 v^r 
82944 


443 h z r* 
1536 
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( 131695 h 2 v 2 r 2 , 38425 h 2 vr 2 , 37681ft 2 ?) 6 ?- _ 400127ft 2 ?) 5 ?- _ 131695ft 2 ?) 4 ?- , 38425ft 2 ?) 3 ?A 

^ 9216 ~ l ~ 36864 5184 165888 55296 ” l " 110592 j axx i-l 


(- 


233 h 2 r 
' 32768 


(- 


(■ 


(- 


1536 

443 h 2 v A 
18432 


505fr 2 i>r 

32768 

505h 2 v 3 _l 
196608 


57161fr 2 i; 7 

995328 


26339 h 2 v 6 
331776 


7685fa 2 ?; £ 

442368 


U r 


27 h 2 
131072 


233 h 2 v 2 i 39 h 2 v 


65536 


262144 j U xx._3 


( 


\ / 2 

24651ft 2 ?) 7 _ 981ft 2 ?) 6 , 3913 u2 6 , 54675ft 2 ?) 5 , 1058589 l 2 2 5 _ 517671ft 2 ?-?) 5 \ , 

v 4096 1024 “*“ 128 n ' u _r 16384 “r 128 U 2048 ^ ^ 

29925 i 1 2 T ,3„,4 , 12879ft 2 ?) 4 , 58695 l2„2, ,4 _ 14715 l2„„4 _ 11659725 /,2„4., 3 , 1764315 u2„3„,3\ , 

v g it I u ~r 1638 4 ~r 128 " ' ft 5!2 "• 1 u 16 u I u ~r 32 '<■ ' u J a xx i _i ~r 

( 54675 ft 2 ?) 3 2588355ft 2 ?- 2 ?) 3 , 273375 h 2 rv 3 , oi i oc l 2„5„,2 89775 ft2 T ,4,,2 , 58695 ft2„3,,2 8019ft 2 ?) 2 A , 

^ 65536 1024 4096 f oiiou/i. /ft 8 a ! u T 32 n I ft 32768 y u-xx^i T 

-^h 2 r 2 v 2 + 3S6 l™l rv2 + 5540535 /l 2 r 6 1 > - + _ 2588355 fc 2 r 3 v + l|6|3ft! 

164025ft 2 ?-?) , ](1 onci,2„6 17955ft 2 ?- 5 , 58695ft 2 ?- 4 14715ft 2 ?- 3 , 6561ft 2 V,, 

32768 4- j -VOVOU I 4 64 128 “*“ 262144 J : 

8019 ft 2 ?- 


(- 


*^xx JL 
1 2 


820125^- 2 r 2 ^ 

4096 


(■ 


( 


(- 


( ! 


(- 


3 8637 ft 2 ? - 2 
4096 


16384 


\ / 2 

24651ft 2 ?' 7 _ 981ft 2 ?) 6 , 3913 l 2 6 __ 54675ft 2 ?) 5 _ 1058589 l 2 2 5 , 517671ft 2 ?-?) 5 A , 

4096 1024 128 ' 1 16384 128 ' v ' 2048 ,1 “**<+! ^ 

29925 , 12879ft 2 ?) 4 , 58695 l2„2, .4 _ 14715 ^2,„„.4 , 11659725 h 2 _4„ 3 _ 1764315 ? 2 „3,,, 3 A ,, , 

v 8 " ' 0 “ft 16384 “ft 128 ^ 512 L u “ft 16 " ' ft 32 " ' ft y u a:a: i+ i ~r 

54675ft 2 ?) 3 , 2588355 ft 2 ?- 2 ?) 3 273375 ft 2 ?-?) 3 , oiiori.2.5,,2 89775 l2„4,, 2 , 58695 l2 t , 3.,2 8019ft 2 ?) 2 A , 

65536 “ft 1024 4096 4 01 100,1 ' 0 s 1 u ‘ 32 n 1 u 32768 ^ “*“ 

44145 ft 2 r2t ,2 + 38637ft 2 r?, 2 _ 55 40535/l 2 r 6 i; + ^f^h 2 r 5 v - ^^h 2 r 4 v + f^h 2 r 3 v - u a:x . + 1 


( 


4096 

820125 h 2 r 2 v , 164025 h 2 rv , 

4096 ' 32768 ' luoyo ^ ' 


(- 


6561^. 2 

262144 


2^2 


38637 h 2 r 
4096 


17955ft 2 r 5 

4 

8019ft 2 r 
16384 


58695ft 2 ?- 4 


64 


14715ftV 

128 


U T 


U x 


200200 ^2^,7 i 63140ft 2 r 6 


( 


( 


\ / ' Z 

1791795^2^6 _ 2313|^ + 63140^^2^ + 654786^2^5 + 2773925^2^3^^ ^ 

188405ft 2 r 4 _ 115675 L2„,2 r 4 _ 818965 ?2„„4 _ 115675 1,2 4_3 _ 818965 l 2„3„3 _ 131695ft 2 r 3 \ 

864 16 n ft ' 576 " ul 48 11 u ' 864 n u ' 13824 ) Uxx i+ 3 

, 188405 i 2 2 3 , 2000635ft 2 ?)?- 3 _ 163793ft 2 ?) 5 r 2 , 188405ft 2 ?) 4 r 2 , 2000635ft 2 ?) 3 r 2 , 443ft 2 ?- 2 

“ft 432 11 y ' “ft 41472 1152 “ft 1728 “ft 82944 “ft 1536 

131695ft 2 ?) 2 r 2 _ 38425ft 2 vr 2 , 37681 ft 2 ?> 6 r , 400127ft 2 ?) 5 r _ 131695ft 2 ?) 4 r _ 38425ft 2 ?) 3 r \ 

9216 36864 “ft 5184 “ft 165888 55296 110592 J 

233ft 2 r I 443ft 2 ?) 2 ?- , 505h 2 vr , 57161ft 2 ?) 7 _ 26339 ft 2 ?) 6 _ 7685 ft 2 ?) 5 \ 

32768 “ft 1536 32768 995328 331776 442368 ^ Ua::E i + | 


( 


(' 


(- 


1536 

443/i 2 f 4 

18432 


505/i 2 t; 3 

196608 


27h 2 

131072 


233h 2 v 2 

65536 

3 „5 


39h 2 v 

262144 




(_ IQOlOftV + H935ft 3 ?- 6 + 155155^3^6 _ 4165ftVi + 11935^3^2^ _ 4^3^) ^ + 

( ■ 65345fr 3 r 4 _ 229075 ^3,,3 r 4 _ 20825 l 3 ? .2 t .4 , 65345/i 3 t;r 4 _ 20825 ^3^4^ _ 22415^ 3 r 3 \ 

6912 288 11 u 1 64 11 u ' ' 1152 192 11 u ' 55296 J u xxx._3 

(i 

(- 


65345ft 3 ?) 3 r 3 


156905ft 3 ?) 3 r 2 
165888 


65345ft 3 ?) 2 r 3 _ 
3456 

22415ft 3 ?) 2 ?- 2 


64 

156905ft 3 ?)?- 

82944 


13069ft 3 ?) 5 ?- 2 

2304 


65345ft 3 ?) 4 ?- 2 


131072 
199ft 3 ?) 5 


995ft 3 ?)?- 2 , 13069ft 3 ?) 6 ?- 31381ft 3 ?) 5 ?- 22415ft 3 ?) 4 ?- \ 

-i- ■- ■— 221184 1 


24576 


41472 


13824 

ft 3 ?) 5 
331776 


199ft 3 r : 
16384 


( 


/ 39ft 3 r I 995ft 3 ?) 3 ?- , 199ft 3 ?) 2 r _ 39 h 3 vr _ 4483ft 3 ?) 7 _ 4483ft 3 ?) 6 \ 

^ 131072 “•“ 73728 16384 65536 1990656 1327104 ^ 

( 


199/i 3 v 


13/i 3 t> 


1048576 131072 



32 


39 h 3 v 2 
262144 


3 h 3 v 
524288 


294912 ' 196608 

135135^ 3 r 6 , 1756755 u3„, r 6 _ 188055 h 3 r 5 , 405405 ^3,,2„5 _ 2068605 ^3.^5 , 156905 h 3 r 4 

4 ' 2 n ui 32 i 4 it u I 16 ri ur 256 


-90090/i 3 r 7 

3447675 ? 3,.3„4 _ 940275 l3 t . 2„4 , 1412145 ? 3 4 _ 313425 ?3 9 .4„3 _ 84015ft 3 r 3 , 470715 l3 ? .3„3 A 

32 11 u ' 64 11 u 1 ' 128 11 u ' 64 n u ‘ 2048 ' 64 11 u 1 J 

■ 156905 u3„,2 r 3 _ 588105 h 3 vr 3 , 282429 u3„,5 r 2 , 156905 l 3,,4„2 , 28917^ 3 r 2 _ 588105/i 3 t; 3 r 2 \ , 

' 128 11 u ' 1024 ' 256 11 u 1 ' 512 11 u 1 ' 16384 2048 J u xxx i _i T 


(- 


(- 
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(- 


252045 h 3 v 2 r 2 
4096 

195 
8192 


144585 h 3 vr 2 
8192 


31381h 3 v 6 r _ 
1536 


(- 


48195ft 3 Y 3 r 


( 4 


117621ft 3 ti 5 r _ 84015ft 3 u 4 r _ 7533ft 3 iA 
4096 8192 131072 J Ux 

5601 h 3 v 7 _ 5601ft 3 Y 6 , 9639 AV\ 

8192 16384 _r 32768 J a xxx._ 


„7 135135ft 3 r 6 

4 


^90090fr 3 r 7 — 4 ~r 2 ,l ul ~ r" 32 

/ 3447675 y3,,3„4 , 940275 l3„.2-.4 , 14121 
^ 32 It u I “r 64 It u I 128 


, 28917ft ! Vr _ 22599ft 3 w _ 

” l " 16384 65536 

539ft 3 t> 4 , 2187ft 3 _ 7533h 3 v 3 _ 7533 h 3 v 2 , 2187ft 3 ?A 
65536 4" 1048576 131072 262144 _r 524288 y axxx i~l 

1756755 y. 3,,,^6 , 188055ft 3 r 5 _ 405405 y 3 ,, 2 _ 5 _ 2068605 l 3„,„5 _ 156905ft 3 r 4 
o ' ” U' i oo 4 it u I it ui 256 


( 


( 4 


(- 


32 4 16 256 ) 1 

1412145 z 1 3 /f . T ,4 , 313425 l 3,,4„3 , 84015fr 3 r 3 , 470715 h 3„,3 r 3\ , 

128 L u ' 64 L u ' 2048 ' 64 L u J u xxx i+ i ' 

t 5~2 156905 y 3 4„2 28917ft 3 r 2 588105ft 3 v 3 r 2 \ , 

' 512 11 u ' 16384 2048 y “*** i+ i 4" 

_ 117621ft 3 Y 5 r I 84015ft 3 Y 4 r , 7533ft 3 r \ , 

4096 8192 131072 y u zxx i+ i f 

48195ft 3 f 3 r _ 28917/t 3 v 2 r _ 22599ft 3 «r _ 5601ft 3 v 7 , 5601ft 3 t> 6 , 9639 ftV 

8192 16384 65536 8192 3 16384 " l " 32768 


156905 y 3 2 3 __ 588105ft 3 «r 3 , 
128 ' 1 1024 


252045ft 3 Y 2 r 2 , 144585ft 3 «r 2 __ 31381ft 3 Y 6 r 

Anna ~~ I 8192 1 Koa 


( 9639 h 3 v 4 _ 2187 h 3 _ 7533 h 3 v 3 , 7533fr 3 u 2 . 

^ 65536 1048576 131072 ' 262144 ' 524288 


2187fc 3 i; 


U'XXX . , 1 
2 


( 


/ 10010ft 3 r 7 _ 11935ft 3 r 6 , 155155 y 3„.„6 , 4165ft 3 r 5 _ 11935 l3,,2„5 _ 45815 u3,„„5\ „ , 

y 3 12 ' is 11 u ' 32 4 u 1 43 11 Ul J Uj xxx i +3 * 


65345 h 3 r 4 _ 229075 l 3 ,,3^,4 , 20825 l 3..2 t .4 
6912 288 ' 64 ^ 


(■ 

( 


65345ft 3 t> 3 r 3 _ 65345ft 3 Y 2 r 3 _ 
1728 3456 


82944 


0 


65345ft 3 t»r 4 , 20825 ft3,,4„3 , 22415ft 3 r 3 \ 

1152 " r 192 “ u 1 ' 55296 ^ 

13069ft 3 f 5 r 2 65345ft 3 v 4 r 2 199 ft 3 r 2 \ 

2304 13824 16384 y u xxx i _ 

13069ft 3 t> 6 r 31381ft 3 -!) 5 r , 22415ft 3 v 4 r\ 

41472 331776 ” l " 221184 ^ a xxx^ 


156905 h 3 v 3 r 2 , 22415ft 3 « 2 r 2 , 995 h 3 vr 2 _ 

165888 ‘ 36864 ' 24576 4I4Y z 331YY0 2: 

995h 3 v 3 r _ 199h 3 v 2 r _ 39 h 3 vr _ 4483 ft 3 t> 7 , 4483ft 3 « 6 \ 
73728 16384 65536 1990656 ” l " 1327104 J 

• 3 v 5 _ 199 h 3 v 4 _ 9ft 3 _ 13ft 3 f 3 , 39 h 3 v 2 j_ 3h 3 v 

Qi? 196608 1048576 131072 4 ” 


39 h 3 r 
131072 


ttxxx ■ 


( 


(- 


( 4 


199ft 3 f 5 _ 
294912 


262144 


4454450ftr 7 

9 


5837755 ftr 6 
36 


90265175 h 6 
54 UUI 


145104575 l 3 4 __ 11173925 y 2 4 , 378708051 
864 uu 1 192 " u 1 ' 20736 

37870805ftt> 2 r 3 _ 13501235ftr 3 _ 110513725ftw 
10368 165888 248832 

13501235ftY» 2 r 2 , 367337 ftr 2 , 714695/tYr 2 
110592 147456 73728 


( 


(- 


U’XXX . 


5837755 y,. 2„5 _ 2234785ftr 5 _ 29020915 y..„5 A „ , 

12 uu 96 144 uu J Ux i- S ' 

I 44422525ft ur 4 _ 11173925 1 4 3 , 44422525fti> 3 r 3 \ 

3456 576 aU 1 ' 5184 J Ux i. 

8884505 hv 5 r 2 , 37870805ftY 4 r 2 _ 110513725fti; 3 r 2 \ 

6912 _l_ 41472 497664 J Ux i-l 

7574161ft-u 6 r _ 22102745ftt> 5 r _ 13501235ftY» 4 r \ 

124416 995328 663552 J 


( 4 


714695^t; 3 r 

221184 


( 


( 4 


367337hv 2 r _ 8105 hr _ 28355 hvr _ 3157535ft» 7 _ 2700247 hv 6 \ 
147456 131072 196608 5971968 3981312 J 


II x. i 
1 2 


142939ft» 5 

884736 


367337 hv 4 
1769472 


28355 ft» 3 _ 8105ft» : 
1179648 262144 


5971968 

1887ft 
1048576 


733hv 

524288 


-1351350/ir 7 + 19438 4 65 ^ r ' 


_ 12818925 to 2 r 4 + 


( 


(- 


1039905/ir 3 _ 27986175^r : 
2048 1024 


( 4 


451251ftr 2 

16384 


_ ^ 695^525^6 + 5831595^2^ _ 2563785^ _ 841^315^5 _ 140280525 ^ 4 ) ^ 

2008335ftr 4 , 60418575; 4 _ 4272975 y,, 4„3 , 20139525 L-.3 r 3 , 2008335 y..2„3 A , 

256 ' 128 lU 64 nu ' 64 ,u 4- l2 8 LV J Ux i-7 

2008335 y 4-2 __ 27986175ftY 3 r 2 _ 3119715ftY 2 r 2 \ , 

512 nu 1 2048 4096 J Ux i-i ” l_ 

5597235ftt> 5 r __ 1039905ftY 4 r , 2939085 hv^r\ , 

4096 8192 4- 81.92 J Ux i - 1 4" 


8817255ftY»r 2 

8192 


12083715 y 5„2 
256 UU 1 

_ 133889 ^ 6 _ 

512 1 


( . 451251ftY 2 r _ 287955 hr _ 2176065ftYr _ 266535ft» 7 _ 69327ft-» 6 , 587817ft» 5 \ 
16384 131072 65536 8192 16384 4- 32768 J 


(l351350/ir 7 - 


( 


' 1 150417ft?; 4 _ 
4" 65536 


725355 hv 3 _ 287955 ftt> 2 , 190269ft , 347733 hv \ 

131072 262144 4" 1048576 4" 524288 I Ux i- 1 


1350 hr 7 - IM3M65ftri + §959^25 hvr e _ 5831595 to 2 r 5 + 25637|5ftri _ 84168315^5 _ 140280525 ^3^ ^ 
12818925^2,, 4 _ 2008335ftr 4 + 604^75^4 + 42^75^4^ + 20139525^3^ _ 2008335^2^^) + 


2563785 ftr 5 _ 84168315 ^ 5 _ 140280525 ^ 3^^ , 

32 16 32 J 


( 4 


1039905ftr 3 _ 
2048 


27986175ftYr 3 , 12083715 . 5 „ 2 _ 2008335 y 4„2 _ 

1024 4- 256 uu 1 512 uu ' 


27986175ft» 3 r 2 

2048 


4096 
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( 451251 hr 2 , 8817255 hvr 2 _ 133889 i 6„ _ 5597235 hv 5 r , 1039905/t» 4 r , 2939085/w 3 r \ , 

^ 16384 8192 512 lv 4096 ^ 8192 ' 8192 J Ux i+± " T " 


(- 


451251/ii> r 
16384 


287955/ir 

131072 


2176065/i t;r 
65536 


266535 hv‘ 


(- 


150417/i 

65536 


725355 hv a 
131072 


8192 

287955 hv 2 190269/i 


69327 hv 6 
16384 


587817 hv 5 
32768 


262144 


1048576 


347733/if \ 
524288 y 


/ 4454450/ir 7 5837755/ir 6 , 90265175 L..-.6 5837755 ^,,2^5 . 2234785/ir 5 29020915 z. 5^ . 

y 9 36 ' 54 ,iUI 12 nu 1 -r 96 144 ,iUI J " r 


145104575 

864 


hv 3 r 4 


11173925 t, 0 ,2^4 


192 


hv z r 4 — 


37870805/ir 

20736 


44422525/ivr 4 


11173925^^3 _l 44422525/iv c V 


2 3 


37870805/it r r' 
10368 

/ , 13501235/wV 


13501235/tr 3 

165888 


2 2 


110592 


367337/ir 

147456 


110513725/i^r 3 

248832 

i 714695/ivr 2 


714695 hv 3 r 
221184 


367337 hvr 
147456 


73728 


8105/ir 

131072 


1 3456 

8884505/1 t; 5 r 2 
6912 

7574161/i t; 6 r 
124416 

28355/i vr 
196608 


576 


_ 37870805/1 t> r 
41472 

22102745fc/i> 5 r , 


5184 


110513725/1^ r 


497664 

13501235/iuV 


995328 

3157535/ro 7 , 
5971968 ' 


2700247hv 6 

3981312 


663552 
U r 


142939 hv 5 367337/iti 4 28355 hv 3 , 8105/i-u 2 


884736 


1769472 


1179648 


262144 


1887/i 

1048576 


733hv 

524288 


M 


X . . 


(74) 


1. cfoS Amplification Factor 

The real part of the c4o3 technique’s amplification factor is shown in Fortran form due to its large size: 


\Re(G)= 

(v* (-588807047879148744*beta + 5029688315970720*beta**3 + 

37059670203005 168640*beta*r - 328971137267483136*beta**3*r - 
1534749197369015130624*beta*r**2 + 
11801835930422937600*beta**3*r**2 + 

4083206537065459 1385600*beta*r**3 - 
5690495 1116416352256*beta**3*r**3 - 
6222948386437471 15008000*beta*r**4 - 
12851851697346192998400*beta**3*r**4 + 
992192170426717962240000*beta*r**5 + 
643479477553473666416640*beta**3*r**5 + 
201463802968439566903541760*beta*r**6 - 
1746 15276987975469 10515200*beta**3*r**6 - 
5286721179261333240859852800*beta*r**7 + 
313381129818851095672258560*beta**3*r**7 + 

71 155 1 140203774977633157 12000*beta*r**8 - 
382230062 1955502779347763200*beta**3*r**8 - 
447028423 148769786829209600000*beta*r**9 + 

28368526 140181720582638796800*beta**3*r**9 - 
1601877496352558815064044339200*beta*r**10 - 
59847788836847349695250432000*beta**3*r**10 + 
588040231785086042650926 12096000*beta*r** 11 - 
1 14558689416 1661438329880576000*beta**3*r** 11 - 
501654460948288553120465682432000*beta*r**12 + 

12480 187407814541778213666816000*beta**3*r** 12 + 
16742359237890354900173 11457280000*beta*r** 13 - 
4442850153354856830 1629603840000*beta**3*r** 13 + 
10856641400139384000*beta*v**2 - 
61032243152201472*beta**3*v**2 - 
86802720302720000 1024*beta*r*v**2 + 
2811635262800179200*beta**3*r*v**2 + 
35316960283429877145600*beta*r**2*v**2 + 
127489092611209003008*beta**3*r**2*v**2 - 
76052585 1586382093352960*beta*r**3*v**2 - 
19737100631473481318400*beta**3*r**3*v**2 + 
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3084002023769 167527936000*beta*r**4*v**2 + 
1039016222710842465976320*beta**3*r**4*v**2 + 
314967737626958823476428800*beta*r**5*v**2 - 
32380786633438227 175833600*beta**3*r**5*v**2 - 
1010746075583 169 1793557094400*beta*r**6*v**2 + 

66 1442776702898338434908 160*beta**3*r**6*v**2 + 
1452547121 1072 1815333253939200*beta*r**7*v**2 - 
8735486 18 185055270 1498163200*beta**3*r**7*v**2 - 
78624890 11 99484763726938 112000*beta*r**8*v**2 + 

65791 903274387010222450278400*beta**3*r**8*v**2 - 
731888710 1806 147288826904576000*beta*r**9*v**2 - 
97419093993513730330066944000*beta**3*r**9*v**2 + 
157738390755256383450315227136000*beta*r**10*v**2 - 
29843923621 18233250213134336000*beta**3*r**10*v**2 - 
959035275737495023494484721664000*beta*r**ll*v**2 + 
23052633926807309847814471680000*beta**3*r**ll*v**2 - 
7925027 1825534920064*beta*v**4 + 
30583554037562880*beta**3*v**4 + 

6421630897 1 17767352320*beta*r*v**4 + 

4749249633700641 1776*beta**3*r*v**4 - 
203859543829363430473728*beta*r**2*v**4 - 
6802413221178936852480*beta**3*r**2*v**4 + 
48620295540803 189 1456000*beta*r**3*v* *4 + 

44235768799 1309237747712*beta**3*r**3*v**4 + 
179634315623296769686241280*beta*r**4*v**4 - 
16792980732045908980531200*beta**3*r**4*v**4 - 
6422236975344074427570585600*beta*r**5*v**4 + 

40097 1483 120935712923320320*beta**3*r**5*v**4 + 
101773594843604206525266001920*beta*r**6*v**4 - 
5885633663530072923753676800*beta**3*r**6*v**4 - 
554704432089795085806 128332800*beta*r**7*v**4 + 
453086068204673985053668 14720*beta**3*r**7*v**4 - 
6059 100162 146023672974095155200*beta*r**8*v**4 - 
26477965 133207704812965068800*beta**3*r**8*v**4 + 
97208 12 19725208786721 18284288000*beta*r**9*v**4 - 
2004030344796760868503368499200*beta**3*r**9*v**4 + 
280563375 107936286720*beta*v**6 + 
3124168881320337408*beta**3*v**6 - 
16539252149112612175872*beta*r*v**6 - 
74160 1712426909368320*beta**3*r*v**6 - 
101728757391537173299200*beta*r**2*v**6 + 

6755703678 1240920834048*beta**3*r**2*v**6 + 
3776794306648631 1089602560*beta*r**3*v**6 - 
3207640886027135012044800*beta**3*r**3*v**6 - 
13961 97232635062331 113472000*beta*r**4*v**6 + 
86683666972156365688012800*beta**3*r**4*v**6 + 

2288 129 17229283902 18239 180800*beta*r**5*v**6 - 
13587 166565037 16645227724800*beta**3*r**5*v**6 - 
125760631653986542616051712000*beta*r**6*v**6 + 
10643200834575207229795860480*beta**3*r**6*v**6 - 
9960 1013449 1996472481520025600*beta*r**7*v**6 - 
1 124625526529472140012421 1200*beta**3*r**7*v**6 - 
358922944231258576896*beta*v**8 - 
23015585822443315200*beta**3*v**8 - 
11941623566508633292800*beta*r*v**8 + 
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3743409734999328817152*beta**3*r*v**8 + 
3210170415656980278804480*beta*r**2*v**8 - 
250908047549 192562278400*beta**3*r**2*v**8 - 
140105518368166183698432000*beta*r**3*v**8 + 
8288905329514443807129600*beta**3*r**3*v**8 + 

22673202341 11262990637465600*beta*r**4*v**8 - 
134105137056035364667392000*beta**3*r**4*v**8 - 
7259350157559910527034982400*beta*r**5*v**8 + 

8054285802 19484723887472640*beta**3*r**5*v**8 - 
495787807221085716480*beta*v**10 + 
69672616691407847424*beta**3*v**10 + 

8860 1438767640 178327552*beta*r*v** 10 - 
7191691599933383639040*beta**3*r*v**10 - 
532 11 696292120 15853568000*beta*r**2*v** 10 + 
316068470842954088448000*beta**3*r**2*v**10 + 

98710 173002797443585146880*beta*r**3*v** 10 - 
5423944737688441690521600*beta**3*r**3*v**10 + 

12839277635627885 19936*bet a* v** 12 - 
79564408586669260800*beta**3*v**12 - 
46335501936495034368000*beta*r*v**12 + 

340898 139 10658 16760320*beta**3*r*v** 12 - 

243*beta* (beta**2* (-57683669631863209448177664000*r**12 + 

1730 15870211 123787302 174720000*r** 13 - 
421906448646144000*r**ll* 

(-21617280715 + 273889751564*v**2) + 

23439247 14700800*r** 10* 

(-399594918611 + 10694645915880*v**2) - 
951268147200*r**8* 

(7496356089035 - 273113326736812*v**2 + 
2137216758218080*v**4) + 

1902536294400*r**9* 

(41595043870467 - 1551530871852500*v**2 + 
7494132128690592*v**4) + 

82575360*r**6* 

(-717986639667879 + 21486595174684016*v**2 - 

161240285877566592* v**4 + 512543000161882368*v**6) - 
660602880*r**7* 

(-1054801734874837 + 32563398611203140*v**2 - 
26410 1846 109232464* v**4 + 583098278454795840*v**6) + 
41287680*r**5* 

(94795433102277 - 3047017849889958*v**2 + 
24033401976772592*v**4 - 78895185920097600*v**6 + 
94232269568938368*v**8) - 

61440*r**4* (3104552759774103 - 108755814337086168*v**2 + 
992484062573122224*v**4 - 3457844668905410560*v**6 + 
5184867573569725440*v**8) - 

24576*r**3* (-271207080146889 + 10250945050778100*v**2 - 
107797789683950562* v**4 + 436899751570237320*v**6 - 
749392863017170400*v**8 + 438028314410897280*v**10) + 
192*r**2* (-829348380467703 + 33079948668167184*v**2 - 
38932083 1684735200* v**4 + 1889296626965601024*v**6 - 
40249845509 10324480*v**8 + 3693650363298488320*v**10) - 
27* (297709173819 - 23981152303080*v**2 + 

342299730399 120*v**4 - 2140251513313536*v**6 + 

682906506 1044480*v* *8 - 10372466907949056*v**10 + 
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6410253269020672*v**12) + 

32*r* (69296470633935 - 3020407834350708*v**2 + 

387259 19708 165280*v**4 - 217803616412366016*v**6 + 
580822312401040896*v**8 - 721684033577183232*v**10 + 
339501476090900480*v**12) ) - 
3* (-775772574024358258 148376576000*r** 12 + 

2438 1327931677930 18 15 1895040000*r** 13 - 
421906448646144000*r**ll* 

(-257254367127 + 3614160898348*v**2) + 

23439247 14700800*r** 10* 

(-3402185718667 + 126271369928760*v**2) + 
66588770304000*r**9* 

(3620684948021 - 401361376717532*v**2 + 
2603560564774752*v**4) - 
951268147200*r**8* 

(-3376274173575 - 1342889495673180*v**2 + 
20112535030790048*v**4) - 
3303014400*r**7* 

(-218271568715125 + 12397919089960596*v**2 - 
2548099260860 13648* v**4 + 1067011446179418432*v**6) + 
82575360*r**6* 

(-2915872934335767 + 52293334784431120*v**2 - 
203709865093847424*v**4 + 2553069348088170240*v**6) + 
41287680*r**5* 

(651620035353411 - 14395064916992790*v**2 + 
447262953551448 16*v**4 - 135736772059360320*v**6 + 
5434355260 19068800* v**8) - 

61440*r**4* (29153895283741815 - 763599752521182360*v**2 + 
4358868307579671984*v**4 - 6185181328054722560*v**6 + 
8606299462809646080*v**8) - 

245760*r**3* (-326013752062737 + 9434326912362198*v**2 - 
73410801524614449*v**4 + 173104319746778292*v**6 - 
767028730 12875760*v**8 + 7731509141205696*v**10) + 
192*r**2* (-12623593377224871 + 384637417525282320*v**2 - 
3498380150575610592*v**4 + 12206778765661843200*v**6 - 
14025471592493326080*v**8 + 4061847949873602560*v**10) \ 

- 27* (6169818240603 - 461835791461800*v**2 + 

5001715839214224*v**4 - 23824855185096960*v**6 + 
57403294649274624*v**8 - 55426485040900096*v**10 + 
66857427434332 16*v** 12) + 

32*r* (1331862912087885 - 44860039946650836*v**2 + 

439976175253469280*v**4 - 1899540445326107328*v**6 + 
3591794938498652160*v**8 - 2682471909002545152*v**10 + 

1 185920327022161920*v**12) ) ) *Cos (beta) - 
972*beta* (beta**2* (-922892772507352674022195200*r**12 + 
21524864780 128553 15595264000*r** 13 - 

84381289729228800*r**ll* (-2365554469 + 24378299660*v**2) + 
23439247 14700800*r** 10* 

(-12158409919 + 243221791632*v**2) - 
761014517760*r**8* 

(306846388688 - 11037490968085*v**2 + 73879224732544*v**4) 
+ 1902536294400*r**9* 

(1555698089817 - 44461683026068*v**2 + 
157343066562528*v**4) + 

82575360*r**6* 
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(-8390260233783 + 400592317432384*v**2 - 
4796853064401408*v**4 + 16887687205248000*v**6) - 
660602880*r**7* 

(-21735341420783 + 918912972202836*v**2 - 
8740846641965616*v**4 + 16939526831464512*v**6) + 
82575360*r**5* 

(319009475466 - 16669199445327*v**2 + 
234442873415960*v**4 - 1144188407488224*v**6 + 
1448247275003712*v**8) - 

6 1440*r**4* (12813557189103 - 703861986264696*v**2 + 
11174990937599664*v**4 - 67997567769399296*v**6 + 

1 52053146 14472 1920* v**8) - 
12288*r**3* (-1500438468915 + 82527016827600*v**2 - 
1422945532585836*v**4 + 10185289876987920*v**6 - 
31139167116948800*v**8 + 28552272358337280*v**10) + 
192*r**2* (-1650077479143 + 89009052247440*v**2 - 
160254766 1864160*v* *4 + 12947818722667776*v**6 - 
4879649656 1053440*v**8 + 76203903203262464*v**10) - 
81* (-59934735 - 9592864344*v**2 + 218226663216*v**4 - 
2010084284160*v**6 + 10031018607360*v**8 - 
26368470 194176*v** 10 + 30903649701888*v**12) + 

64*r* (46120044303 - 2849042564514*v**2 + 

533451 16022976*v**4 - 459076826220768*v**6 + 
2011462884761472*v**8 - 4254001425607680*v**10 + 
2909977010382848*v**12) ) - 
3* (-20133583894660502098280448000*r**12 + 

514232990 15953439407472640000*r** 13 - 
421906448646144000*r**ll* 

(-9679278267 + 101865765716*v**2) + 

23439247 14700800*r** 10* 

(-245765287109 + 4905031931328*v**2) - 
1902536294400*r**8* 

(2959552682415 - 96219577428774*v**2 + 

604433839 1981 12*v**4) + 

1902536294400*r**9* 

(33223356231995 - 906542266954940*v**2 + 
3243439678639776*v**4) - 
660602880*r**7* 

(-620539208944357 + 22640836761551100*v**2 - 

18880962 1496456976* v**4 + 341796645237347520*v**6) + 
82575360*r**6* 

(-288543121430805 + 11671815463754336*v**2 - 

116857257267352320*v**4 + 364001677713148416*v**6) + 
82575360*r**5* 

(13104388796445 - 584971692371685*v**2 + 
6766961878348840*v**4 - 27939797535271200*v**6 + 
31561275189465024*v**8) - 

61440*r**4* (612130122784485 - 29493473801105160*v**2 + 
389745499 167568656* v**4 - 1958865473739317248*v**6 + 
36976070 15636997120*v**8) - 
12288*r**3* (-80726104925487 + 4027934222870400*v**2 - 
593607229 13632644*v**4 + 352703853912613680*v**6 - 
887189420132478400*v**8 + 681034036144062720*v**10) + 
192*r**2* (-98476779983445 + 4922372676072816*v**2 - 
77715621186994080*v**4 + 534442122526535424*v**6 - 
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1675966626530069760*v**8 + 2177600498496274432*v**10) - 
27* (14928433569 - 2207820187512*v**2 + 

40585421655216*v**4 - 328743522743040*v**6 + 
1427223068762880*v**8 - 3169036058834944*v**10 + 
30556327098 16320*v** 12) + 

128*r* (1651783345704 - 88594168691835*v**2 + 

1464721895747088*v**4 - 11067394312967760*v**6 + 
41628975096003264*v**8 - 73850069684021760*v**10 + 
43397862838787072*v**12) ) ) *Cos (2*beta) - 
3*beta* (-6565413079503 + 398408459938326331671969792000*r**12 - 
90608808 110 1032576086 179840000*r** 13 + 110821561318728*v**2 - 
401211 159351 120*v**4 - 4508563972688640*v**6 + 

586565046 16268544*v**8 - 257715918703392768*v**10 + 
442123540826525696*v**12 + 

421906448646144000*r**ll* 

(-199054845767 + 2182021247692*v**2) - 
2343924714700800*r**10* 

(-4729245686083 + 103656419886264*v**2) + 

951268147200*r**8* 

(69844313562865 - 2972871503160924*v**2 + 

2294662 1498684576*v**4) - 
1902536294400*r**9* 

(531609461763455 - 17239248316624180*v**2 + 

66 1342440 14058528*v**4) - 

82575360*r**6* (-1385166227803471 + 87325784946246928*v**2 - 
135 156578086 161 1392*v* *4 + 5758348552837801728*v**6) + 
660602880*r**7* (-4865642730162041 + 256056788770319940*v**2 - 
297 1496980680499728*v**4 + 6653161396905291840*v**6) - 
41287680*r**5* (70947893367747 - 5326699377817030*v**2 + 
103227072160039792*v**4 - 641970137232605760*v**6 + 
979355319246328704*v**8) + 

430080*r**4* (115071039657225 - 10705019299468200*v**2 + 
257430127321849680*v**4 - 2138251218314289152*v**6 + 
6115711802512051200*v**8) + 

49152*r**3* (-7713649675629 + 1227142504767870*v**2 - 
38765752221783573*v**4 + 414931844050079140*v**6 - 
17326088052 18292400*v* *8 + 2050361230337542080*v**10) - 
192*r**2* (36136825502337 + 1550210230422288*v**2 - 

101250428823753696*v**4 + 1447856822418607872*v**6 - 
80 127784737808 11520*v* *8 + 16755880172081340416*v**10) - 
32*r* (-11208197155347 + 166109466520332*v**2 + 

2472521938752096*v**4 - 61285385287929024*v**6 + 
446619625746905088*v**8 - 1360377446820857856*v**10 + 

1 18654622759 1983104*v** 12) + 

3*beta**2* (-659 10850909 125 1205072486400*r** 12 + 

15233205 18370749755424768000*r** 13 - 

84381289729228800*r**ll* (-1612681903 + 18059678780*v**2) + 
2343924714700800*r**10* (-7484112451 + 168099152424*v**2) + 
1902536294400*r**9* 

(821151171723 - 27290457944596*v**2 + 

10728 1230487200*v**4) - 
190253629440*r**8* 

(526084952279 - 22934452411900*v**2 + 
181725049039072*v**4) + 

82575360*r**6* 
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(-1956547557303 + 127748838546160*v**2 - 
2032 17603671 1552*v* *4 + 8885105774274816*v**6) - 
660602880*r**7* 

(-7130914417565 + 385061401604292*v**2 - 
4585461138962640*v**4 + 10544008313745984*v**6) + 
41287680*r**5* 

(93887635293 - 7526238878022*v**2 + 

15102 127 1377264*v**4 - 963450698392896*v**6 + 
1510243631355264*v**8) - 

61440*r**4* (892544804055 - 99481400362008*v**2 + 
2549241380372400*v**4 - 21869195086548992*v**6 + 
64186134339855360*v**8) - 
24576*r**3* (-1149062193 + 2780159930100*v**2 - 
104119336285650*v**4 + 1175524010390280*v**6 - 
5058029287 141600*v* *8 + 6153750656565120*v**10) + 
192*r**2* (110404472985 + 77801877648*v**2 - 

1 19366959569 120*v**4 + 1954509147636480*v**6 - 
11342305865137920*v**8 + 24358626486591488*v**10) + 
3*(3711406797 - 67873440024*v**2 + 358447305456*v**4 + 

1 1 17039880448* v**6 - 24498416113920*v**8 + 
117600703752192*v**10 - 209771938164736*v**12) + 

32*r* (-21509011593 + 420604058700*v**2 + 

129 1629672864* v**4 - 74872746083520*v**6 + 
605030639961600*v**8 - 1917925378615296*v**10 + 
1704777524703232*v**12) ) ) *Cos (3*beta) - 
729* (bet a**2* (105741 172866477573305008 128000*r** 12 - 
345588754177326433582448640000*r**13 + 
4640970935107584000*r**ll* 

(-2850552659 + 44081501412*v**2) - 
16407473002905600*r**10* 

(-37497573689 + 2176321827240*v**2) + 

951268147200*r**8* 

(-8243164839537 + 23669349438860*v**2 + 
1861865235800416*v**4) - 
9512681472000*r**9* 

(-4347617889887 - 244038396953428*v**2 + 
2266640436651360*v**4) + 

82575360*r**6* 

(201824348917275 + 5393607906635120*v**2 - 

1 17193940856745344*v**4 + 40426634189395200*v**6) + 
3303014400*r**7* 

(114111407047925 - 3788820566188332*v**2 + 
5949837097675920*v**4 + 93817478273572416*v**6) - 
41287680*r**5* 

(94093784713071 - 1250846573090382*v**2 - 
8969352529977392*v**4 + 49587992868653760*v**6 + 

985389 1643771520*v* *8) - 

61440*r**4* (-4931360757517671 + 111162807540305640*v**2 - 
34148800967512 1904* v**4 - 1349570909297766400*v**6 + 
333968255797360 1280*v**8) - 

122880*r**3* (115311651991713 - 3205582644247668*v**2 + 

21 18680583 1200 198*v**4 - 22636216114826184*v**6 - 
82502 104862326240*v**8 + 85388234245087872*v**10) + 
192*r**2* (2172688760787939 - 68755345498995120*v**2 + 
594437324079273696*v**4 - 1732248770187859200*v**6 + 
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657607468274369280*v**8 + 1979673708999024640*v**10) - 
27* (-1009995530787 + 74309862642840*v**2 - 

875005671090384*v**4 + 4293338196259584*v**6 - 
9746141498802432*v**8 + 7380223941482496*v**10 + 

197158 16598 11840*v** 12) + 

32*r* (-213827634433485 + 7809056613423684*v**2 - 

78990312626485344*v**4 + 322917377690563776*v**6 - 
50435086 1926325760* v**8 + 150948224433948672*v**10 + 
13505346093854720*v**12) ) + 

2* (-579164488159548445053419520000*r**12 + 
1961191715319373579576934400000*r**13 - 
4219064486461440000*r**ll* 

(-15359554547 + 260703519934*v**2) + 
11719623573504000*r**10* 

(-83020302737 + 14746647361080*v**2) - 
9512681472000*r**8* 

(-9761508015143 + 129985977196304*v**2 + 
572297552058848*v**4) + 

38050725888000*r**9* 

(-16875641979912 - 157709025385835*v**2 + 
2841462403951248*v**4) - 
6606028800*r**7* 

(992837565004693 - 29231979456070830*v**2 + 

131 19 1705522068768* v**4 + 125738952524375520*v**6) - 
412876800*r**6* 

(-538711671355491 + 30466191272258128*v**2 - 
325903046690708736* v**4 + 482733945503753472*v**6) - 
41287680*r**5* 

(-103170194733375 - 8146923188365390*v**2 + 

1 93505373764547392* v**4 - 736520632101652800*v**6 + 

3 1039973178322 1760* v**8) + 

307200*r**4* (-3133463798505879 + 29163429763252848*v**2 + 
637275042406935408* v**4 - 5656716299520763904*v**6 + 
98081 11688628126720*v**8) + 

61440*r**3* (948011830571583 - 19617646470228720*v**2 + 
43577030 10 1694544*v**4 + 681912210566310480*v**6 - 
2834559992969425280*v**8 + 2110158314000513280*v**10) - 
960*r**2* (2158755861170673 - 53898431669866944*v**2 + 
346419096157008864*v**4 - 173053798755784704*v**6 - 
3763428220024482560*v**8 + 6893215787690860544*v**10) + 
9* (-20129692202145 + 1449671858802000*v**2 - 

1229608688 1752592*v**4 + 47340492487649280*v**6 - 
72658653626085 120*v**8 - 40528416929239040*v**10 + 

1 65266499427307520* v** 12) - 
160*r* (-278158668457275 + 7481171537805564*v**2 - 
59885380704290976*v**4 + 176103003953732544*v**6 - 
42506522 17341696*v* *8 - 619379876900419584*v**10 + 
363292515836551 168* v** 12) ) ) *Sin(beta) + 

1458* (931491792905523037376348 16000*r** 12 - 
246877865688655000523243520000*r**13 + 

421906448646144000*r**ll* (-43287796421 + 462914161996*v**2) - 
11719623573504000*r**10* 

(-216435349655 + 4359492372864*v**2) - 
104639496 192000*r**9* 

(2698620485907 - 72933563676988*v**2 + 264833023123872*v**4) \ 
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+ 1902536294400*r**8* 

(13805511831569 - 434288957577530*v**2 + 

2680609668773344* v**4) + 

33030 14400*r**7* 

(-613864546602751 + 21323646693005844*v**2 - 

169854088665737520*v**4 + 301270771112410176*v**6) - 
412876800*r**6* (-302987318486055 + 11629586218331552*v**2 - 
10957440 1523831040*v**4 + 327268567196702208*v**6) - 
82575360*r**5* (72551512578735 - 3088646615010355*v**2 + 
33597781 165796696*v**4 - 131158083360499680*v**6 + 
14264882365 10 19840*v**8) + 

61440*r**4* (3532323834269955 - 163994687251347960*v**2 + 
2053511071558037616*v**4 - 9719048608637419520*v**6 + 
17333519957424491520*v**8) + 

61440*r**3* (-95983707035637 + 4666317548673600*v**2 - 
65918071771476876*v**4 + 370758980823898896*v**6 - 
8769 1348841 17 14880*v**8 + 633480120919505664*v**10) - 
960*r**2* (-119780938536327 + 5879215360055376*v**2 - 
89877487849922016*v**4 + 591656152031882496*v**6 - 
1756770906870959872*v**8 + 2154121759149359104*v**10) - 
640*r* (207 1 112379030 - 108494532089289*v**2 + 

1746481645544688*v**4 - 12777458758845552*v**6 + 
46084793355740736*v**8 - 77624201846410752*v**10 + 
435485888 11277312*v** 12) + 

9*(350489177685 - 41944107421080*v**2 + 

745001267346288*v**4 - 5869666090817280*v**6 + 
24618475207745280*v**8 - 52183701319014400*v**10 + 
47443699365826560*v**12) + 
beta**2* (-5851842503476545885044736000*r**12 + 

1458451269 1569717902573568000*r** 13 - 
84381289729228800*r**ll* 

(-14478492031 + 149040099780*v**2) + 

23439247 14700800*r** 10* 

(-75253922725 + 1470778685712*v**2) - 
3805072588800*r**8* 

(448677732615 - 14739871368323*v**2 + 92102778883808*v**4) 
+ 1902536294400*r**9* 

(10238612744539 - 276767674825180*v**2 + 
967663477977504*v**4) - 
660602880*r**7* 

(-181466996077885 + 6823946603640540*v**2 - 
579329 14257676944*v**4 + 104446520276894400*v**6) + 
82575360*r**6* 

(-80535955555197 + 3395324450176640*v**2 - 
35318602303764480*v**4 + 111698694999201792*v**6) + 
82575360*r**5* 

(3488621369112 - 162446990875029*v**2 + 

1973806767302536*v**4 - 8439247978353312*v**6 + 
965569875 1530432*v* *8) - 

430080*r**4* (22370894453907 - 1116354114342936*v**2 + 
15490999 132628976*v**4 - 81679616621400064*v**6 + 
159694719017303040*v**8) - 

12288*r**3* (-20043453482109 + 1025042979865680*v**2 - 
15744856565515140*v**4 + 98344407722704560*v**6 - 
25975297 1261963200* v**8 + 206870355639978240*v**10) + 
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192*r**2* (-23881268931597 + 1213649688770352*v**2 - 
19805244850217376*v**4 + 142185814593488640*v**6 - 
468574408032318720*v**8 + 637535391567757312*v**10) - 
27* (2583570897 - 510424962264*v**2 + 9739448713008*v**4 - 
81264819767040*v**6 + 365709873537792*v**8 - 
849889633572864*v**10 + 863173582524416*v**12) + 

64*r* (775381807791 - 42530703638166*v**2 + 

723105192345984*v**4 - 5648368335055776*v**6 + 
22146137322791040*v**8 - 41206726443156480*v**10 + 
25169881674348544*v**12) ) ) *Sin(2*beta) + 

(9*beta**2* (258734708055 - 15549235181138349398163456000*r**12 + 
3562307143874688134519 1936000*r** 13 - 4550593234104*v**2 + 
20409795982224*v**4 + 127896243169536*v**6 - 
20 12012824953600* v**8 + 9192060007925760*v**10 - 
16060874463268864*v**12 - 
928194187021516800*r**ll* 

(-3497258887 + 38725811540*v**2) + 

23439247 14700800*r** 10* 

(-180757714199 + 4008278601624*v**2) + 
13317754060800*r**9* 

(2868920831899 - 94143364395740*v**2 + 
365360529862944*v**4) - 
951268147200*r**8* 

(2606643331919 - 112245625081044*v**2 + 

8774062502941 12*v**4) + 

82575360*r**6* 

(-50147232969059 + 3213399676226960*v**2 - 
50398337480389248*v**4 + 217442907736898304*v**6) - 
660602880*r**7* 

(-179206790012575 + 9548202735783300*v**2 - 

112209081625691568*v**4 + 254474776093496640*v**6) + 
41287680*r**5* 

(2494786113855 - 192837447947102*v**2 + 

37987 18875072080*v**4 - 23918871981926208*v**6 + 

3697 124366907 1488*v**8) - 

61440*r**4* (26320552610031 - 2637942453019944*v**2 + 
65270894920980336*v**4 - 550491961326555136*v**6 + 
1594432797576299520*v**8) - 
24576*r**3* (-319858286049 + 80848429244100*v**2 - 
2742568464175470*v**4 + 30075460967614120*v**6 - 
1 27397218441488800* v**8 + 152803760846808960*v**10) + 
192*r**2* (1994151317589 + 31633066726704*v**2 - 
3388133391429600*v**4 + 51329238364212480*v**6 - 
290300 1635786 12480* v**8 + 614902888963170304*v**10) + 
32*r* (-470940516531 + 8160191550684*v**2 + 

638619537 14784*v**4 - 2080336736786112*v**6 + 
15855413162964480*v**8 - 49194214246434816*v**10 + 
43308340658917376*v**12) ) - 

2*00166464489525 - 1842098813469762639567519744000*r**12 + 
4164992816402544456923873280000*r**13 - 
487728092905920*v**2 + 1308808254098448*v**4 + 
26450188396715520*v**6 - 303662123843800320*v**8 + 
1292370351419 105280*v** 10 - 2180785886053765120*v**12 - 
5062877383753728000*r**ll* 

(-77330254127 + 840198168667*v**2) + 
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35 1588707205 12000*r** 10* 

(-1484489830757 + 32201219462232*v**2) - 
7610145177600*r**8* 

(42023175324339 - 1769435026490695*v**2 + 

1350041 1156203664*v**4) + 

19025362944000*r**9* 

(253006881664231 - 8116017220136826*v**2 + 
30818687380278336*v**4) - 
1 189085 18400*r**7* 

(-131705920213382 + 6852511340086811*v**2 - 
78586799399 187768*v**4 + 173909698940190224*v**6) + 
412876800*r**6* 

(-1370490234992227 + 85171108096236432*v**2 - 

1301996090061407232*v**4 + 5482204198604107008*v**6) + 
41287680*r**5* 

(359545292830305 - 26358525922696670*v**2 + 
503356757273686944* v**4 - 3094050621272811840*v**6 + 
466319189941 1308800*v**8) - 

552960*r**4* (479017085154345 - 42174611799274080*v**2 + 
99037 17242 10410416*v**4 - 8113035789389404160*v**6 + 
2293029 1222525 137920*v**8) - 
61440*r**3* (-41456131042059 + 5228930118173760*v**2 - 
1565268401864931 12* v**4 + 1641186719654641744*v**6 - 
6763031918331279360*v**8 + 7902816330307804416*v**10) + 
960*r**2* (20048562602793 + 2098046110342752*v**2 - 

106703534850323616*v**4 + 1459078409555891712*v**6 - 
7925337652 125707008*v**8 + 16375351377067941888*v**10) \ 

+ 1440*r* (-1068273665775 + 12808707849756*v**2 + 
338133889156320*v**4 - 7104681800882496*v**6 + 

49965 166436758272*v**8 - 149745393445836800*v**10 + 
129453832 195702784* v** 12) ) ) *Sin(3*beta) ) )/ 

2 . 19 1220846 1633946e 17 


The imaginary part of the amplification factor for the c4o3 method is shown here in Fortran form due to 
its length: 


\Im(G)= 

(v*(-588807047879148744*beta + 5029688315970720*beta**3 + 

37059670203005 168640*beta*r - 328971137267483136*beta**3*r - 
1534749197369015130624*beta*r**2 + 
11801835930422937600*beta**3*r**2 + 

4083206537065459 1385600*beta*r**3 - 
5690495 1116416352256*beta**3*r**3 - 
622294838643747115008000*beta*r**4 - 
12851851697346192998400*beta**3*r**4 + 
992192170426717962240000*beta*r**5 + 
643479477553473666416640*beta**3*r**5 + 
201463802968439566903541760*beta*r**6 - 
1746 15276987975469 10515200*beta**3*r**6 - 
5286721179261333240859852800*beta*r**7 + 
313381129818851095672258560*beta**3*r**7 + 

71 155 1 140203774977633157 12000*beta*r**8 - 
382230062 1955502779347763200*beta**3*r**8 - 
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447028423 148769786829209600000*beta*r**9 + 

28368526 140181720582638796800*beta**3*r**9 - 
160 18774963525588 15064044339200*beta*r** 10 - 
59847788836847349695250432000*beta**3*r**10 + 
588040231785086042650926 12096000*beta*r** 11 - 
1145586894161661438329880576000*beta**3*r**ll - 
501654460948288553120465682432000*beta*r**12 + 

12480 1874078145417782 13666816000*beta**3*r** 12 + 
16742359237890354900173 11457280000*beta*r** 13 - 
4442850153354856830 1629603840000*beta**3*r** 13 + 
10856641400139384000*beta*v**2 - 
61032243152201472*beta**3*v**2 - 
86802720302720000 1024*beta*r*v**2 + 
2811635262800179200*beta**3*r*v**2 + 
35316960283429877145600*beta*r**2*v**2 + 
127489092611209003008*beta**3*r**2*v**2 - 
76052585 1586382093352960*beta*r**3*v**2 - 
19737100631473481318400*beta**3*r**3*v**2 + 
3084002023769 167527936000*beta*r**4*v**2 + 
1039016222710842465976320*beta**3*r**4*v**2 + 
314967737626958823476428800*beta*r**5*v**2 - 
32380786633438227 175833600*beta**3*r**5*v**2 - 
1010746075583169 1793557094400*beta*r**6*v**2 + 

66 1442776702898338434908 160*beta**3*r**6*v**2 + 
1452547121 1072 1815333253939200*beta*r**7*v**2 - 
8735486 18 185055270 1498163200*beta**3*r**7*v**2 - 
78624890 11 99484763726938 112000*beta*r**8*v**2 + 

65791 903274387010222450278400*beta**3*r**8*v**2 - 
731888710 1806 147288826904576000*beta*r**9*v**2 - 
97419093993513730330066944000*beta**3*r**9*v**2 + 
157738390755256383450315227136000*beta*r**10*v**2 - 
2984392362 118233250213134336000*beta**3*r**10*v**2 - 
959035275737495023494484721664000*beta*r**ll*v**2 + 
23052633926807309847814471680000*beta**3*r**ll*v**2 - 
79250271825534920064*beta*v**4 + 
30583554037562880*beta**3*v**4 + 

642 1630897 117767352320*beta*r*v**4 + 

4749249633700641 1776*beta**3*r*v**4 - 
203859543829363430473728*beta*r**2*v**4 - 
6802413221178936852480*beta**3*r**2*v**4 + 
48620295540803 189 1456000*beta*r**3*v**4 + 

44235768799 1309237747712*beta**3*r**3*v**4 + 
179634315623296769686241280*beta*r**4*v**4 - 
16792980732045908980531200*beta**3*r**4*v**4 - 
6422236975344074427570585600*beta*r**5*v**4 + 

40097 1483 120935712923320320*beta**3*r**5*v**4 + 
101773594843604206525266001920*beta*r**6*v**4 - 
5885633663530072923753676800*beta**3*r**6*v**4 - 
554704432089795085806 128332800*beta*r**7*v**4 + 
453086068204673985053668 14720*beta**3*r**7*v**4 - 
6059 100162 146023672974095155200*beta*r**8*v**4 - 
26477965 133207704812965068800*beta**3*r**8*v**4 + 
97208 12 19725208786721 18284288000*beta*r**9*v**4 - 
2004030344796760868503368499200*beta**3*r**9*v**4 + 
280563375 107936286720*beta*v**6 + 
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3124168881320337408*beta**3*v**6 - 
16539252149112612175872*beta*r*v**6 - 
74160 1712426909368320*beta**3*r*v**6 - 
101728757391537173299200*beta*r**2*v**6 + 

6755703678 1240920834048*beta**3*r**2*v**6 + 

3776794306648631 1089602560*beta*r**3*v**6 - 
3207640886027135012044800*beta**3*r**3*v**6 - 
13961 97232635062331 113472000*beta*r**4*v**6 + 
86683666972156365688012800*beta**3*r**4*v**6 + 

2288 129 17229283902 18239 180800*beta*r**5*v**6 - 
13587 166565037 16645227724800*beta**3*r**5*v**6 - 
125760631653986542616051712000*beta*r**6*v**6 + 
10643200834575207229795860480*beta**3*r**6*v**6 - 
9960 1013449 199647248 1520025600*beta*r**7*v**6 - 
1 124625526529472140012421 1200*beta**3*r**7*v**6 - 
358922944231258576896*beta*v**8 - 
23015585822443315200*beta**3*v**8 - 
11941623566508633292800*beta*r*v**8 + 

37434097349993288 17 152*beta**3*r*v**8 + 
3210170415656980278804480*beta*r**2*v**8 - 
250908047549 192562278400*beta**3*r**2*v**8 - 
140105518368166183698432000*beta*r**3*v**8 + 
8288905329514443807129600*beta**3*r**3*v**8 + 

22673202341 11262990637465600*beta*r**4*v**8 - 
134105137056035364667392000*beta**3*r**4*v**8 - 
7259350157559910527034982400*beta*r**5*v**8 + 

8054285802 19484723887472640*beta**3*r**5*v**8 - 
495787807221085716480*beta*v**10 + 
69672616691407847424*beta**3*v**10 + 

8860 1438767640 178327552*beta*r*v** 10 - 
7191691599933383639040*beta**3*r*v**10 - 
532 11 696292120 15853568000*beta*r**2*v** 10 + 
316068470842954088448000*beta**3*r**2*v**10 + 

98710 173002797443585146880*beta*r**3*v** 10 - 
5423944737688441690521600*beta**3*r**3*v**10 + 
12839277635627885 19936*bet a* v** 12 - 
79564408586669260800*beta**3*v**12 - 
46335501936495034368000*beta*r*v**12 + 

340898 139 10658 16760320*beta**3*r*v** 12 - 
243*beta*(beta**2*(-57683669631863209448177664000*r**12 + 

1730 15870211 123787302 174720000*r** 13 - 
421906448646144000*r**ll* 

(-21617280715 + 273889751564*v**2) + 

23439247 14700800*r** 10* 

(-399594918611 + 10694645915880*v**2) - 
951268147200*r**8* 

(7496356089035 - 273113326736812*v**2 + 
2137216758218080*v**4) + 

1902536294400*r**9* 

(41595043870467 - 1551530871852500*v**2 + 
7494132128690592*v**4) + 

82575360*r**6* 

(-717986639667879 + 21486595174684016*v**2 - 

161240285877566592* v**4 + 512543000161882368*v**6) - 
660602880*r**7* 
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(-1054801734874837 + 32563398611203140*v**2 - 
26410 1846 109232464* v**4 + 583098278454795840*v**6) + 
41287680*r**5* 

(94795433102277 - 3047017849889958*v**2 + 
24033401976772592*v**4 - 78895185920097600*v**6 + 
94232269568938368*v**8) - 

61440*r**4* (3104552759774103 - 108755814337086168*v**2 + 
992484062573122224*v**4 - 3457844668905410560*v**6 + 
5184867573569725440*v**8) - 

24576*r**3* (-271207080146889 + 10250945050778100*v**2 - 
107797789683950562*v**4 + 436899751570237320*v**6 - 
749392863017170400*v**8 + 438028314410897280*v**10) + 
192*r**2* (-829348380467703 + 33079948668167184*v**2 - 
389320831684735200*v**4 + 1889296626965601024*v**6 - 
40249845509 10324480*v**8 + 3693650363298488320*v**10) - 
27* (297709173819 - 23981152303080*v**2 + 

342299730399 120*v**4 - 2140251513313536*v**6 + 

682906506 1044480*v* *8 - 10372466907949056*v**10 + 
6410253269020672*v**12) + 

32*r* (69296470633935 - 3020407834350708*v**2 + 

387259 19708 165280*v**4 - 217803616412366016*v**6 + 
580822312401040896*v**8 - 721684033577183232*v**10 + 
339501476090900480*v**12) ) - 
3* (-775772574024358258 148376576000*r** 12 + 

2438 1327931677930 18 15 1895040000*r** 13 - 
421906448646144000*r**ll* 

(-257254367127 + 3614160898348*v**2) + 

23439247 14700800*r** 10* 

(-3402185718667 + 126271369928760*v**2) + 
66588770304000*r**9* 

(3620684948021 - 401361376717532*v**2 + 
2603560564774752*v**4) - 
951268147200*r**8* 

(-3376274173575 - 1342889495673180*v**2 + 
20112535030790048*v**4) - 
3303014400*r**7* 

(-218271568715125 + 12397919089960596*v**2 - 
2548099260860 13648* v**4 + 1067011446179418432*v**6) + 
82575360*r**6* 

(-2915872934335767 + 52293334784431120*v**2 - 
203709865093847424*v**4 + 2553069348088170240*v**6) + 
41287680*r**5* 

(651620035353411 - 14395064916992790*v**2 + 
447262953551448 16*v**4 - 135736772059360320*v**6 + 
5434355260 19068800* v**8) - 

61440*r**4* (29153895283741815 - 763599752521182360*v**2 + 
4358868307579671984*v**4 - 6185181328054722560*v**6 + 
8606299462809646080*v**8) - 

245760*r**3* (-326013752062737 + 9434326912362198*v**2 - 
73410801524614449*v**4 + 173104319746778292*v**6 - 
767028730 12875760*v**8 + 7731509141205696*v**10) + 
192*r**2* (-12623593377224871 + 384637417525282320*v**2 - 
3498380150575610592*v**4 + 12206778765661843200*v**6 - 
14025471592493326080*v**8 + 4061847949873602560*v**10) \ 

- 27* (6169818240603 - 461835791461800*v**2 + 
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5001715839214224*v**4 - 23824855185096960*v**6 + 
57403294649274624*v**8 - 55426485040900096*v**10 + 
66857427434332 16*v** 12) + 

32*r* (1331862912087885 - 44860039946650836*v**2 + 

439976175253469280*v**4 - 1899540445326107328*v**6 + 
3591794938498652160*v**8 - 2682471909002545152*v**10 + 

1 185920327022161920*v**12) ) ) *Cos (beta) - 
972*beta* (beta**2* (-922892772507352674022195200*r**12 + 
21524864780 128553 15595264000*r** 13 - 

84381289729228800*r**ll* (-2365554469 + 24378299660*v**2) + 
23439247 14700800*r** 10* 

(-12158409919 + 243221791632*v**2) - 
761014517760*r**8* 

(306846388688 - 11037490968085*v**2 + 73879224732544*v**4) 
+ 1902536294400*r**9* 

(1555698089817 - 44461683026068*v**2 + 
157343066562528*v**4) + 

82575360*r**6* 

(-8390260233783 + 400592317432384*v**2 - 
4796853064401408*v**4 + 16887687205248000*v**6) - 
660602880*r**7* 

(-21735341420783 + 918912972202836*v**2 - 
8740846641965616*v**4 + 16939526831464512*v**6) + 
82575360*r**5* 

(319009475466 - 16669199445327*v**2 + 

234442873415960*v**4 - 1144188407488224*v**6 + 
1448247275003712*v**8) - 

61440*r**4* (12813557189103 - 703861986264696*v**2 + 
11174990937599664*v**4 - 67997567769399296*v**6 + 

1 52053146 14472 1920* v**8) - 
12288*r**3* (-1500438468915 + 82527016827600*v**2 - 
1422945532585836*v**4 + 10185289876987920*v**6 - 
31139167116948800*v**8 + 28552272358337280*v**10) + 
192*r**2* (-1650077479143 + 89009052247440*v**2 - 
160254766 1864160*v**4 + 12947818722667776*v**6 - 
4879649656 1053440*v**8 + 76203903203262464*v**10) - 
81* (-59934735 - 9592864344*v**2 + 218226663216*v**4 - 
2010084284160*v**6 + 10031018607360*v**8 - 
26368470 194176*v** 10 + 30903649701888*v**12) + 

64*r* (46120044303 - 2849042564514*v**2 + 

533451 16022976*v**4 - 459076826220768*v**6 + 
2011462884761472*v**8 - 4254001425607680*v**10 + 
2909977010382848*v**12) ) - 
3* (-20133583894660502098280448000*r**12 + 
51423299015953439407472640000*r**13 - 
421906448646144000*r**ll* 

(-9679278267 + 101865765716*v**2) + 

23439247 14700800*r** 10* 

(-245765287109 + 4905031931328*v**2) - 
1902536294400*r**8* 

(2959552682415 - 96219577428774*v**2 + 

604433839 1981 12*v**4) + 

1902536294400*r**9* 

(33223356231995 - 906542266954940*v**2 + 
3243439678639776*v**4) - 
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660602880*r**7* 

(-620539208944357 + 22640836761551100*v**2 - 

18880962 1496456976* v**4 + 341796645237347520*v**6) + 
82575360*r**6* 

(-288543121430805 + 11671815463754336*v**2 - 

116857257267352320*v**4 + 364001677713148416*v**6) + 
82575360*r**5* 

(13104388796445 - 584971692371685*v**2 + 
6766961878348840*v**4 - 27939797535271200*v**6 + 
31561275189465024*v**8) - 

61440*r**4* (612130122784485 - 29493473801105160*v**2 + 
389745499 167568656* v**4 - 1958865473739317248*v**6 + 
36976070 15636997120*v**8) - 
12288*r**3* (-80726104925487 + 4027934222870400*v**2 - 
593607229 13632644*v**4 + 352703853912613680*v**6 - 
887189420132478400*v**8 + 681034036144062720*v**10) + 
192*r**2* (-98476779983445 + 4922372676072816*v**2 - 
77715621186994080*v**4 + 534442122526535424*v**6 - 
1675966626530069760*v**8 + 2177600498496274432*v**10) 
27* (14928433569 - 2207820187512*v**2 + 

40585421655216*v**4 - 328743522743040*v**6 + 
1427223068762880*v**8 - 3169036058834944*v**10 + 
30556327098 16320*v** 12) + 

128*r* (1651783345704 - 88594168691835*v**2 + 

1464721895747088*v**4 - 11067394312967760*v**6 + 
41628975096003264*v**8 - 73850069684021760*v**10 + 
43397862838787072*v**12) ) ) *Cos (2*beta) - 
3*beta* (-6565413079503 + 398408459938326331671969792000*r**12 
90608808 110 1032576086 179840000*r** 13 + 110821561318728*v**2 
401211159351120*v**4 - 4508563972688640*v**6 + 

586565046 16268544*v**8 - 257715918703392768*v**10 + 
442123540826525696*v**12 + 

421906448646144000*r**ll* 

(-199054845767 + 2182021247692*v**2) - 
2343924714700800*r**10* 

(-4729245686083 + 103656419886264*v**2) + 
951268147200*r**8* 

(69844313562865 - 2972871503160924*v**2 + 

2294662 1498684576*v**4) - 
1902536294400*r**9* 

(531609461763455 - 17239248316624180*v**2 + 

66 1342440 14058528*v**4) - 

82575360*r**6* (-1385166227803471 + 87325784946246928*v**2 - 
135 156578086 161 1392*v**4 + 5758348552837801728*v**6) + 
660602880*r**7* (-4865642730162041 + 256056788770319940*v**2 
2971496980680499728*v**4 + 6653161396905291840*v**6) - 
41287680*r**5* (70947893367747 - 5326699377817030*v**2 + 
103227072160039792*v**4 - 641970137232605760*v**6 + 
979355319246328704*v**8) + 

430080*r**4* (115071039657225 - 10705019299468200*v**2 + 
257430 12732 1849680*v**4 - 2138251218314289152*v**6 + 
6115711802512051200*v**8) + 

49152*r**3* (-7713649675629 + 1227142504767870*v**2 - 
38765752221783573*v**4 + 414931844050079140*v**6 - 
17326088052 18292400*v* *8 + 2050361230337542080*v**10) - 


NASA/TM— 2005-213976 


45 



192*r**2* (36136825502337 + 1550210230422288*v**2 - 

101250428823753696*v**4 + 1447856822418607872*v**6 - 
80 127784737808 11520*v* *8 + 16755880172081340416*v**10) - 
32*r* (-11208197155347 + 166109466520332*v**2 + 

247252 1938752096* v**4 - 61285385287929024*v**6 + 
446619625746905088*v**8 - 1360377446820857856*v**10 + 

1 18654622759 1983104*v** 12) + 

3*beta**2* (-659 10850909 125 1205072486400*r** 12 + 

15233205 18370749755424768000*r** 13 - 

84381289729228800*r**ll* (-1612681903 + 18059678780*v**2) + 
2343924714700800*r**10* (-7484112451 + 168099152424*v**2) + 
1902536294400*r**9* 

(821151171723 - 27290457944596*v**2 + 

10728 1230487200*v**4) - 
190253629440*r**8* 

(526084952279 - 22934452411900*v**2 + 
181725049039072*v**4) + 

82575360*r**6* 

(-1956547557303 + 127748838546160*v**2 - 
2032 17603671 1552*v**4 + 8885105774274816*v**6) - 
660602880*r**7* 

(-7130914417565 + 385061401604292*v**2 - 
4585461 138962640*v**4 + 10544008313745984*v**6) + 
41287680*r**5* 

(93887635293 - 7526238878022*v**2 + 

15102 127 1377264*v**4 - 963450698392896*v**6 + 
1510243631355264*v**8) - 

61440*r**4* (892544804055 - 99481400362008*v**2 + 
2549241380372400*v**4 - 21869195086548992*v**6 + 
64186134339855360*v**8) - 
24576*r**3* (-1149062193 + 2780159930100*v**2 - 
104119336285650*v**4 + 1175524010390280*v**6 - 
5058029287 141600*v* *8 + 6153750656565120*v**10) + 
192*r**2* (110404472985 + 77801877648*v**2 - 

1 19366959569 120*v**4 + 1954509147636480*v**6 - 
11342305865137920*v**8 + 24358626486591488*v**10) + 
3*(3711406797 - 67873440024*v**2 + 358447305456*v**4 + 

1 1 17039880448* v**6 - 24498416113920*v**8 + 

1 17600703752 192*v** 10 - 209771938164736*v**12) + 

32*r* (-21509011593 + 420604058700*v**2 + 

129 1629672864* v**4 - 74872746083520*v**6 + 
605030639961600*v**8 - 1917925378615296*v**10 + 
1704777524703232*v**12) ) ) *Cos (3*beta) - 
729* (bet a**2* (105741 172866477573305008 128000*r** 12 - 
345588754177326433582448640000*r**13 + 
4640970935107584000*r**ll* 

(-2850552659 + 44081501412*v**2) - 
16407473002905600*r**10* 

(-37497573689 + 2176321827240*v**2) + 

951268147200*r**8* 

(-8243164839537 + 23669349438860*v**2 + 
1861865235800416*v**4) - 
9512681472000*r**9* 

(-4347617889887 - 244038396953428*v**2 + 
2266640436651360*v**4) + 
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82575360*r**6* 

(201824348917275 + 5393607906635120*v**2 - 

1 17193940856745344*v**4 + 40426634189395200*v**6) + 
3303014400*r**7* 

(114111407047925 - 3788820566188332*v**2 + 
5949837097675920*v**4 + 93817478273572416*v**6) - 
41287680*r**5* 

(94093784713071 - 1250846573090382*v**2 - 
8969352529977392*v**4 + 49587992868653760*v**6 + 

985389 1643771520*v* *8) - 

61440*r**4* (-4931360757517671 + 111162807540305640*v**2 - 
341488009675121904*v**4 - 1349570909297766400*v**6 + 
333968255797360 1280*v**8) - 

122880*r**3* (115311651991713 - 3205582644247668*v**2 + 

21 18680583 1200 198*v**4 - 22636216114826184*v**6 - 
82502 104862326240*v**8 + 85388234245087872*v**10) + 
192*r**2* (2172688760787939 - 68755345498995120*v**2 + 
594437324079273696*v**4 - 1732248770187859200*v**6 + 
657607468274369280*v**8 + 1979673708999024640*v**10) - 
27* (-1009995530787 + 74309862642840*v**2 - 

875005671090384*v**4 + 4293338196259584*v**6 - 
9746141498802432*v**8 + 7380223941482496*v**10 + 

197158 16598 11840*v** 12) + 

32*r* (-213827634433485 + 7809056613423684*v**2 - 

78990312626485344*v**4 + 322917377690563776*v**6 - 
50435086 1926325760* v**8 + 150948224433948672*v**10 + 
13505346093854720*v**12) ) + 

2* (-579164488159548445053419520000*r**12 + 
1961191715319373579576934400000*r**13 - 
4219064486461440000*r**ll* 

(-15359554547 + 260703519934*v**2) + 
11719623573504000*r**10* 

(-83020302737 + 14746647361080*v**2) - 
9512681472000*r**8* 

(-9761508015143 + 129985977196304*v**2 + 
572297552058848*v**4) + 

38050725888000*r**9* 

(-16875641979912 - 157709025385835*v**2 + 
2841462403951248*v**4) - 
6606028800*r**7* 

(992837565004693 - 29231979456070830*v**2 + 

131 19 1705522068768* v**4 + 125738952524375520*v**6) - 
412876800*r**6* 

(-538711671355491 + 30466191272258128*v**2 - 
325903046690708736*v**4 + 482733945503753472*v**6) - 
41287680*r**5* 

(-103170194733375 - 8146923188365390*v**2 + 

1 93505373764547392* v**4 - 736520632101652800*v**6 + 

3 1039973178322 1760* v**8) + 

307200*r**4* (-3133463798505879 + 29163429763252848*v**2 + 
637275042406935408* v**4 - 5656716299520763904*v**6 + 
98081 11688628126720*v**8) + 

61440*r**3* (948011830571583 - 19617646470228720*v**2 + 
43577030 10 1694544*v**4 + 681912210566310480*v**6 - 
2834559992969425280*v**8 + 2110158314000513280*v**10) - 
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960*r**2* (2158755861170673 - 53898431669866944*v**2 + 
346419096157008864*v**4 - 173053798755784704*v**6 - 
3763428220024482560*v**8 + 6893215787690860544*v**10) + 
9* (-20129692202145 + 1449671858802000*v**2 - 

1 229608688 1752592*v**4 + 47340492487649280*v**6 - 
72658653626085 120*v**8 - 40528416929239040*v**10 + 

1 65266499427307520* v** 12) - 
160*r* (-278158668457275 + 7481171537805564*v**2 - 
59885380704290976*v**4 + 176103003953732544*v**6 - 
42506522 17341696*v* *8 - 619379876900419584*v**10 + 
363292515836551 168* v** 12) ) ) *Sin(beta) + 

1458* (93149 1792905523037376348 16000*r** 12 - 
246877865688655000523243520000*r**13 + 

421906448646144000*r**ll* (-43287796421 + 462914161996*v**2) - 
11719623573504000*r**10* 

(-216435349655 + 4359492372864*v**2) - 
104639496 192000*r**9* 

(2698620485907 - 72933563676988*v**2 + 264833023123872*v**4)\ 

+ 1902536294400*r**8* 

(13805511831569 - 434288957577530*v**2 + 

2680609668773344* v**4) + 

33030 14400*r**7* 

(-613864546602751 + 21323646693005844*v**2 - 

169854088665737520*v**4 + 301270771112410176*v**6) - 
412876800*r**6* (-302987318486055 + 11629586218331552*v**2 - 
10957440 1523831040*v**4 + 327268567196702208*v**6) - 
82575360*r**5* (72551512578735 - 3088646615010355*v**2 + 
33597781 165796696*v**4 - 131158083360499680*v**6 + 
14264882365 10 19840*v**8) + 

61440*r**4* (3532323834269955 - 163994687251347960*v**2 + 
2053511071558037616*v**4 - 9719048608637419520*v**6 + 
17333519957424491520*v**8) + 

61440*r**3* (-95983707035637 + 4666317548673600*v**2 - 
65918071771476876*v**4 + 370758980823898896*v**6 - 
8769134884117 14880*v**8 + 633480120919505664*v**10) - 
960*r**2* (-119780938536327 + 5879215360055376*v**2 - 
89877487849922016*v**4 + 591656152031882496*v**6 - 
1756770906870959872*v**8 + 2154121759149359104*v**10) - 
640*r* (207 1 112379030 - 108494532089289*v**2 + 

1746481645544688*v**4 - 12777458758845552*v**6 + 
46084793355740736*v**8 - 77624201846410752*v**10 + 

435485888 11277312*v** 12) + 

9* (350489177685 - 41944107421080*v**2 + 

745001267346288*v**4 - 5869666090817280*v**6 + 
24618475207745280*v**8 - 52183701319014400*v**10 + 
47443699365826560*v**12) + 
beta**2* (-5851842503476545885044736000*r**12 + 

1458451269 1569717902573568000*r** 13 - 
84381289729228800*r**ll* 

(-14478492031 + 149040099780*v**2) + 

23439247 14700800*r** 10* 

(-75253922725 + 1470778685712*v**2) - 
3805072588800*r**8* 

(448677732615 - 14739871368323*v**2 + 92102778883808*v**4) 
+ 1902536294400*r**9* 
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(10238612744539 - 276767674825180*v**2 + 
967663477977504*v**4) - 
660602880*r**7* 

(-181466996077885 + 6823946603640540*v**2 - 
579329 14257676944*v**4 + 104446520276894400*v**6) + 
82575360*r**6* 

(-80535955555197 + 3395324450176640*v**2 - 
35318602303764480*v**4 + 111698694999201792*v**6) + 
82575360*r**5* 

(3488621369112 - 162446990875029*v**2 + 

1973806767302536*v**4 - 8439247978353312*v**6 + 
965569875 1530432*v* *8) - 

430080*r**4* (22370894453907 - 1116354114342936*v**2 + 
15490999 132628976*v**4 - 81679616621400064*v**6 + 

1 596947190 17303040* v**8) - 

12288*r**3* (-20043453482109 + 1025042979865680*v**2 - 
1 5744856565515 140*v**4 + 98344407722704560*v**6 - 
25975297 1261963200* v**8 + 206870355639978240*v**10) + 
192*r**2* (-23881268931597 + 1213649688770352*v**2 - 
19805244850217376*v**4 + 142185814593488640*v**6 - 
468574408032318720*v**8 + 637535391567757312*v**10) - 
27* (2583570897 - 510424962264*v**2 + 9739448713008*v**4 - 
81264819767040*v**6 + 365709873537792*v**8 - 
849889633572864*v**10 + 863173582524416*v**12) + 

64*r* (775381807791 - 42530703638166*v**2 + 

723105192345984*v**4 - 5648368335055776*v**6 + 
22146137322791040*v**8 - 41206726443156480*v**10 + 
25169881674348544*v**12) ) ) *Sin(2*beta) + 

(9*beta**2* (258734708055 - 15549235181138349398163456000*r**12 + 
3562307143874688134519 1936000*r** 13 - 4550593234104*v**2 + 
20409795982224*v**4 + 127896243169536*v**6 - 
20 12012824953600* v**8 + 9192060007925760*v**10 - 
16060874463268864*v**12 - 
928194187021516800*r**ll* 

(-3497258887 + 38725811540*v**2) + 

23439247 14700800*r** 10* 

(-180757714199 + 4008278601624*v**2) + 
13317754060800*r**9* 

(2868920831899 - 94143364395740*v**2 + 
365360529862944*v**4) - 
951268147200*r**8* 

(2606643331919 - 112245625081044*v**2 + 

8774062502941 12*v**4) + 

82575360*r**6* 

(-50147232969059 + 3213399676226960*v**2 - 
50398337480389248*v**4 + 217442907736898304*v**6) - 
660602880*r**7* 

(-179206790012575 + 9548202735783300*v**2 - 

112209081625691568*v**4 + 254474776093496640*v**6) + 
41287680*r**5* 

(2494786113855 - 192837447947102*v**2 + 

37987 18875072080*v**4 - 23918871981926208*v**6 + 

3697 124366907 1488*v**8) - 

61440*r**4* (26320552610031 - 2637942453019944*v**2 + 
65270894920980336*v**4 - 550491961326555136*v**6 + 
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1594432797576299520*v**8) - 
24576*r**3* (-319858286049 + 80848429244100*v**2 - 
2742568464175470*v**4 + 30075460967614120*v**6 - 
1273972 18441488800* v**8 + 152803760846808960*v**10) + 
192*r**2* (1994151317589 + 31633066726704*v**2 - 
3388133391429600*v**4 + 51329238364212480*v**6 - 
290300 1635786 12480* v**8 + 614902888963170304*v**10) + 
32*r* (-470940516531 + 8160191550684*v**2 + 

63861953714784*v**4 - 2080336736786112*v**6 + 
15855413162964480*v**8 - 49194214246434816*v**10 + 
43308340658917376*v**12) ) - 

2*(30166464489525 - 1842098813469762639567519744000*r**12 + 
4164992816402544456923873280000*r**13 - 
487728092905920*v**2 + 1308808254098448*v**4 + 
26450188396715520*v**6 - 303662123843800320*v**8 + 
1292370351419 105280*v** 10 - 2180785886053765120*v**12 - 
5062877383753728000*r**ll* 

(-77330254127 + 840198168667*v**2) + 

35 1588707205 12000*r** 10* 

(-1484489830757 + 32201219462232*v**2) - 
7610145177600*r**8* 

(42023175324339 - 1769435026490695*v**2 + 

1350041 1156203664*v**4) + 

19025362944000*r**9* 

(253006881664231 - 8116017220136826*v**2 + 
30818687380278336*v**4) - 
1 189085 18400*r**7* 

(-131705920213382 + 6852511340086811*v**2 - 
78586799399187768*v**4 + 173909698940190224*v**6) + 
412876800*r**6* 

(-1370490234992227 + 85171108096236432*v**2 - 

1301996090061407232*v**4 + 5482204198604107008*v**6) + 
41287680*r**5* 

(359545292830305 - 26358525922696670*v**2 + 
503356757273686944* v**4 - 3094050621272811840*v**6 + 
466319189941 1308800*v**8) - 

552960*r**4* (479017085154345 - 42174611799274080*v**2 + 
99037 17242 10410416* v**4 - 8113035789389404160*v**6 + 
2293029 1222525 137920*v**8) - 
61440*r**3* (-41456131042059 + 5228930118173760*v**2 - 
1565268401864931 12* v**4 + 1641186719654641744*v**6 - 
6763031918331279360*v**8 + 7902816330307804416*v**10) + 
960*r**2* (20048562602793 + 2098046110342752*v**2 - 

106703534850323616* v**4 + 1459078409555891712*v**6 - 
7925337652 125707008*v**8 + 16375351377067941888*v**10) \ 

+ 1440*r* (-1068273665775 + 12808707849756*v**2 + 
338133889156320*v**4 - 7104681800882496*v**6 + 

49965 166436758272*v**8 - 149745393445836800*v**10 + 
129453832 195702784* v** 12) ) ) *Sin(3*beta) ) )/ 

2 . 19 1220846 1633946e 17 
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